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Abstract 

In order to understand the nonlinear stability of many types of time-periodic travelling waves on un- 
bounded domains, one must overcome two main difficulties: the presence of embedded neutral eigenvalues 
and the time-dependence of the associated linear operator. This problem is studied in the context of time- 
periodic Lax shocks in systems of viscous conservation laws. Using spatial dynamics and a decomposition 
into separate Floquet eigenmodes, it is shown that the linear evolution for the time-dependent operator can 
be represented using a contour integral similar to that of the standard time-independent case. By decompos- 
ing the resulting Green's distribution, the leading order behavior associated with the embedded eigenvalues 
is extracted. Sharp pointwise bounds are then obtained, which are used to prove that time-periodic Lax 
shocks are linearly and nonlinearly stable under the necessary conditions of spectral stability and minimal 
multiplicity of the translational eigenvalues. The latter conditions hold, for example, for small-oscillation 
time-periodic waves that emerge through a supercritical Hopf bifurcation from a family of time-independent 
Lax shocks of possibly large amplitude. 
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1 Introduction 



We are interested in solutions to partial differential equations on unbounded domains that are time-periodic 
in an appropriate comoving frame and spatially asymptotic to constants. Such solutions, which we refer to as 
time-periodic waves, are therefore of the form 

u{x,t) = u{x — ct,t) with lim u{y,t)=u±, u{y,t + T)=u{y,t), (1-1) 

y — ^±00 

where c is the speed of the wave, and T > is the time period. We refer to travelling waves u{x — ci) as time- 
independent waves because they become stationary in the comoving frame. Time-periodic solutions to partial 
differential equations arise in many applications. Examples are pulsating, spinning, and cellular instabilities in 
detonation waves, which have been investigated in [21 [H], and interfaces that connect wave trains in reaction- 
diffusion systems, which have been analyzed, for instance, in |23| . The purpose of this paper is to prove that 
spectral stability of time-periodic waves in viscous conservation laws implies their nonlinear stability. In the 
examples given above, two difficulties arise that prevent us from applying standard theory: the linearization 
about the wave is time-periodic, and its Floquet spectrum contains essential spectrum up to the origin. 

We will focus our analysis on time-periodic Lax shocks, which arise in systems of viscous conservation laws. This 
setting is useful not only because of its relation to real physical models that exhibit periodic phenomena, but also 
because the underlying structure within these systems provides intuition and aids in their analysis. As discussed 
in [24l [26l EH EH] , time-periodic viscous Lax shocks may arise through Hopf bifurcations from stationary shock 
waves. Spectral stability of the bifurcating waves has been treated in [53], where it was shown that they are 
spectrally stable if the bifurcation is supercritical and spectrally unstable if the bifurcation is subcritical. Our 
goal is to prove nonlinear stability of arbitrary time-periodic Lax shocks, possibly far away from any stationary 
shocks. In the stationary case, pointwise Green's function estimates have proved to be very useful in establishing 
nonlinear stability; see, for instance, [71[T3]. We shall use ideas from spatial dynamics, and in particular the 
exponential-dichotomy theory for ill-posed elliptic problems developed in [TT] EU E3] > to extend this approach to 
the time-periodic case. 

We now state our hypotheses and results in detail. Consider a time-periodic viscous shock profile t) = 
u{x — ct,t), as in ( |1.1[ ), of a parabolic system of conservation laws 

ut + fiu)^ = u^^, x£R, mgM^, (1.2) 

where f £ . Working in a coordinate system that moves along with the shock and appropriately rescaling 
time, we may, without loss of generality, consider a standing profile u{x^t) with minimal temporal period 27r. 
Accordingly, we take c—Q and T = 2tt from now on. We shall assume that the profile u{x, t) is a Lax p-shock: 

Hypothesis (HI) The ordered eigenvalues af < . . . < of fu{u±) are real, distinct, and nonzero, and there 
is a number p £{!,... , N} so that a^_p < < ajj_p^i and ajf_p^i < < aj^_p^2- 

The eigenvalues a J of Ju{u±) determine the characteristics x = a^t oi the linear system 

vt + fu{u±)vx = 0, 

whose general solution is a linear combination of solutions of the form v^x — a^t). We say that the characteristics 
corresponding to a^T < and a^t > are outgoing, while the characteristics associated with a^T > and < 
are called incoming, and refer the reader to Figure [I] for an illustration. For future use, we denote by l^ and 
the left and right eigenvectors of fu{u±) associated with the characteristic speeds and normalize these 
vectors so that {l~,r'^) = (^^,r^) = Sjk- 
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Figure 1: The N -\-\ incoming and N — 1 outgoing far-field characteristics x = a^t of the linearization about a 
Lax p-shock for a system of N viscous conservation laws are illustrated for p = 2 and TV = 3. 



Definition 1.1 Throughout this paper, we will often use the notation a^^^ and a^^^ to denote the outgoing 
characteristics < and > o/ u_ and respectively. Similarly, and a^^ denote the incoming 
characteristics a^T > and a'j^ < 0. We often write 

E-E E-E E-E E-E 

to emphasize the interpretation of these sums in terms of characteristics. Similar notation will also be used for 
the associated left and right eigenvectors. 

In order to state our spectral assumptions, we linearize ( |1.2| about the time-periodic shock u{x,t) and obtain 

ut = u^x - [fu{u{x,t))u]r,, (1.3) 

where fu{u{x,t)) is 27r-periodic in t. Because of the time-periodicity of the linear operator on the right-hand 
side, the appropriate characterization of spectral stability is in terms of its Floquet spectrum. Let <f>27r denote 
the time-27r map associated with the semiflow of (1.3l posed on L^(]R,M^). The Floquet spectrum is defined as 

E = {o- e C : c^""" e spectrum of <^2^}- (1-4) 

It follows from Hypothesis |(H1)] that Ux{x,t) and ut{x,t) decay exponentially to zero as \x\ —> oo, uniformly 
in time; see [24: . Thus, the origin a = lies in the spectrum E and belonging to it are two exponentially 
decaying eigenfunctions, namely Ux and ut, which refiect the translation-invariance of the underlying equations 
with respect to space and time. We will assume that the origin is the only part of the spectrum that lies in the 
closed right half-plane other than the spectrum at iZ that results from the fact that Floquet exponents are not 
unique. 

Definition 1.2 The profile u is said to he spectrally stable if 

(SI) The time-2TT map ^2it associated with {1.3) on L^(M,]R^) has no eigenvalues in {Re cr > 0} \ iZ. 



(52) Equation (1.3) has exactly two linearly independent localized solutions with period 2tt in t. 

(53) The outgoing characteristics and the jump [u] := — 7i_ across the shock u are linearly independent: 

det(rf , . . . , r]^_p, [u], r+_^2-p' ..■,rjj)^0. (1.5) 



(S4) The formal adjoint equation 

wt = Wxx + fl{u{x, t))Wx 

has a non-constant, hounded, time 2'K-periodic solution w{x,t). Setting tp2{x,t) 
that 

{'ilj2{x,t),ut{x,t)) dtdx ^ 



w{x,~t), we assume 
(1.6) 



and limx^±oo 4^2{x,t) =: ^± € L^, where denote the suhspaces o/M^ spanned hy the left eigenvectors 
l^ of fu{u±) associated with the incoming characteristics a^. 
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Hypothesis (S2) guarantees that the geometric multiphcity of the zero eigenvalue is minimal. Condition (S3) 
is the same as for stationary Lax shocks, where it is known as the Liu-Majda condition; see [33] for more 
details. Furthermore, (S3) implies that there is a unique, up to scalar multiples, nonzero ipi G M.^ so that ipi 
is perpendicular to the outgoing characteristic eigenvectors that appear in the determinant (1.5 1. This constant 
V'l is the adjoint eigenfunction associated with Ux, and (S3) implies that 

r27T I 

(^i,?2,(a;,t))da;di= (V^i,[?2]) ^0. (1.7) 

I 

Given (S2), conditions (S3)-(S4) furthermore imply that zero is a root of order two of an appropriate Evans 
function. The integrals in (1.6) and (1.7 1 can be interpreted as the Melnikov integrals associated with speed and 



frequency for the spatial-dynamics existence problem for u{x,t); see [22l |24] . 

Our first result, paralleling the theory of the time-independent case fl41l33j. is the equivalence of linearized and 
spectral stability. 



Theorem 1 Assume (HI) and pick p>0, then spectral stability is equivalent to linearized stability in C] 
(that is, each solution of (1.3) with initial data in C\ converges in this space to spanj-ita;, Wt} as t —> ooj. 

Our next theorem is the main result of this paper: it asserts that spectral stability implies nonlinear stability 
and gives detailed pointwise estimates for how perturbations decay as i — > cxd. Before we can state this result, 
we need additional notation. Let 



^ gauss 



{x,t) 



E 



-\x-a^^^t\^ /Mt 



E 



-|£c-a+„,t|^/Aft 



,{x,t) 



E 



E 



V^+¥\+t ^l + \x-a-^,t\ ^ ^l + |x-a+ti| 



{l + \x-a-^t\ + ^/ty2 {l + \x-a+t\ + Vi)l 
where M > is a sufficiently large constant, and define 



X ^ [a'^t.a'l^t] 



to be the characteristic function of the characteristic cone [a^ t, a^t]. 



Theorem 2 Define the weighted norm \\v\\fj3 ■= |1(1 -t- a:^)iu||/f3 . // (HI) and spectral stability (S1)-(S4) hold, 
then the profile u is nonlinearly stable with respect to initial perturbations vq for which Hwolli/J, is sufficiently 
small. More precisely, there exist constants C > and S > such that, for each vq with \\vq\\h3 < S, there exist 
functions {q,T){t) and constants (q^,T^) so that, for all x CzR and t>0, we have 

\u{x,t) - u{x - q^ - q{t),t - n - T{t))\ < C||^;o||^r3 [6'gauss + X^'innor + (1 - x)^outGr] (a;, i) (1-8) 
\\u{-,t) - u{- ~ q^ ~ q{t),t - - T{t))\\H3 < C\\vq\\h3 



\{q.,n)\ + (1 + tp\{q,T){t)\ + (1 + t)\{q,T){t)\ < C\\vo\\m, 



where u{x,t) denotes the solution of (1.2) with initial data uq(x) = u{x,0) + vo{x). 



The pointwise bound (1.8) yields as a corollary the sharp decay rate 

\\u{-,t) - u{- - - q{t),t - r, - T{t))\\LP < C||?;o||h3 (1 + t)-^'-^--p' 



1 < p < oo. 
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The statement of Theorem[2]can be understood as follows. First, for initial data sufficiently close to the underlying 
time-periodic wave, solutions to the full nonlinear equation (1.2 1 will converge to an appropriate space and time 
translate of the wave. In addition, the perturbation to the solution profile will be of a certain form, given by 
the functions ^gauss, ^'inncr, ^ud 6'outer, as it decays to zero. The function flgauss consists of Gaussians that move 
along the outgoing characteristics at speeds aj < and aj" > 0. The function f?inncr accounts for the nonlinear 
interactions that occur within the characteristic cone [a^<, aj^t] encoded in the characteristic function x- The 
algebraically decaying tail of the initial data, outside the characteristic cone, is captured by the function Pouter- 

Theorem [2] is obtained as a consequence of detailed estimates of the solution operator of the linearization 

ut = u^x - [fu{u{x,t))u]rc. (1.9) 



More precisely, we shall derive pointwise bounds of the Green's distribution Q(x,t\y, s) of (1.9 1, which is its 
fundamental solution, defined as the solution at (x, t) of ( 1.9 1 with initial data 5{x — y) at time s, where 5{x — y) 
is the delta-function centered at y. Recall that u{x, t) converges to the constant vectors u± exponentially fast as 
X ±cx), and that a^^^ and aj^ are, respectively, the outgoing and incoming characteristics associated with u±. 
We denote by errfn(z) := ^ e~^ d^ the error function. 

Theorem 3 Under the assumptions of Theorem^ the Green's distribution Q{x,t;y,s) associated with the lin- 
earized system ( 1.9 I can be written as Q = £i + £2 + G so that the following is true: There are positive constants 
Tj, C , and M so that, for y < and t > s, we have 

£i(x,t;y, s) = Ux{x,t)iii{y,s,t) (1-10) 
£2{x,t;y,s) = Ut{x,t)Tr2{y, s,t) (1.12) 

for appropriate functions Zj,in(y, s) that are 2TT-periodic in s and satisfy 

mi,,{y,s)\ + \d:i2,Ay,s)\ < C, \dfd^lT,Ay,^)\ + \d:d^l2,,^{y,s)\ < Ce-"!^! (1.14) 
for < |a| < 1 and I < (3 < 2, and 

\d^g{x,t-y,s)\ < Ce-''(l^-^l+(*-^» 

k e-(^-y-<^~ (t-o)f /Mit-s)^-vx+ 



+C ((t - s)-^ + lale-"!^!) ^(t - s) 

\ a~ 

_L V V (f ^ s)-h-^-^^-''^Mt~'>)-\v/<ut\)f/M{t-s) -r,x+ 



-i^-<utiit-'>)-\y/<^J))/M{t-s)^-rix- 



for < a < 1. In the above, x^ :— max{0,±a;} > denotes the positive/negative part of x, and Xj{^) ^s the 
indicator function of the interval J. Symmetric bounds hold for y > and t > s. Furthermore, if we replace 
u{x, t) with u{x — qi,,t — r*), the estimates above remain true uniformly for (g*, r*) in any compact set. 

The terms £1 and £2 correspond to the translational eigenmodes in space and time, respectively. The error 
functions in ( |1.11[ ) and (1.131 record the effects of information that gets transported along the incoming char- 
acteristics, which determine the ultimate space and time translate of the wave to which the perturbed solution 
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converges. The term Q encapsulates higher order terms arising from the parts of the spectrum to the left of the 
imaginary axis, including the continuous spectrum near zero, and from higher order terms associated with the 
translational eigenmodes. Overall, the description of the Green's distribution in Theorem [3] is exactly analogous 
to that in [7] for the time-independent case, with the addition of the new term £2- Thus, in Theoremjs] we have 
effectively performed a "time-asymptotic conjugation" to the time-independent case, analogous to the usual Flo- 
quet conjugation for finite-dimensional differential equations, but carried out only on the low-frequency modes 
important for time-asymptotic behavior. 

We remark that we have made no effort to state minimal regularity assumptions; instead, we chose to maintain an 
argument structure that generalizes easily to the case of quasilinear and partially parabolic viscosity {B(u)ux)x 
as treated in [T51 [55] . In the constant- viscosity case treated here, could be replaced by in Theorem [2] 
and the short-time theory of §6.2| and §6.3| which is based on energy estimates, could be replaced by simpler L°° 
theory based on integral equations and Picard iteration as in [7, §5]. We expect that our results will extend to 
the quasilinear, partially parabolic case as in the closely related analyses of E5] . 

We have focused here on pure Lax shocks, which are the most common type of viscous shocks, and the only 
type arising in standard gas dynamics. Other types of stationary viscous shocks include undercompressive, 
overcompressive, and mixed- type profiles: pure over- and undercompressive profiles arise in magnetohydrody- 
namics (MHD) and phase-transitional models, but mixed under-overcompressive profiles are also possible, as 
described in 133] . The nonlinear stability of stationary non-Lax viscous shocks has recently been addressed 
in [7]. One key difference is that, in the pure Lax and the undercompressive case, the neutral eigenvalues re- 
sult only from space translates of the profile u, whereas in other cases they also involve deformations of m; see 
[HI [T31 15U1 ISn 1311 133] for further discussion. Our results for the case of time-periodic viscous shocks carry over 
to nonclassical over-, under-, and mixed over-undercompressive waves, with virtually no changes in the analysis, 
provided we substitute for the stability condition (S2) the condition that there exist I eigenmodes at cr = 0, 
where £ is the dimension of the manifold of time-periodic travelling-wave connections between u_ and it+ and 
for (S3)-(S4) the more fundamental condition 



det 



2ir 



{i>j{x,t), cj)k{x,t)) dx d< 



for bases ipj and (/f^ of left and right genuine eigenmodes of tr = 0, where j,k — !,...,£. With this change, the 
analysis goes through essentially unchanged; see fZ] for the time-independent case. 

The plan for the remainder of the paper is as follows. In fj2] we give a detailed explanation of the method by which 
we will prove the main results of the paper, focusing on the intuition behind the ideas, rather than the technical 
details. Section [3] contains the linear theory we develop for equation (1.9 1 using spatial dynamics, including the 



contour integral representation of the linear evolution, which provides the framework for the decomposition of the 
Green's distribution in ^|4] Using this decomposition. Theorems [T| and [3] are proved in f|5] while f|6] contains the 
nonlinear analysis and the proof of Theorem |2] We end this paper in S|7] with a brief summary and a discussion 
of open problems. 



2 Outline of the method 

We now give a detailed outline of the method used to obtain the results in ^ focusing on the intuition and main 
ideas involved. Rigorous justification can be found in Sj3]-Sj6] below. 
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2.1 Nonlinear stability of stationary shocks 



We first recall the techniques that can be used to study the nonlinear stability of stationary Lax shocks u(x) of 



(1.2 1. Taking a solution of the form u(x,t) — u(x) + v{x,t), the resulting equation for the perturbation v is 

vt ^ Cv + Q{v,Vx)x, v{-,0) = vo, 

where 

J^V = - [fu{u{x))v]x, Q{V, Vx) = ~f{u + v) + f{u) + fu{u)v. 

One way to understand the linearized equation 

Vt — Cv 

is to take its Laplace transform and obtain 

Xi) — vq ~ Lv. 

This equation can be solved using the resolvent operator, and taking the inverse Laplace transform of the result 
leads to the standard representation 



e*^ 



^ / e*^[A-/:]-MA 



27ri 

of the linear semigroup, where F is a curve in the complex plane that lies entirely in the resolvent set of H,. 
For viscous shocks, the essential spectrum of £ contains the origin. Hence, we cannot easily derive good decay 
estimates for the linear semigroup, because we cannot move the contour F into the origin, where analyticity 
of the resolvent [A — breaks down. Instead, we exploit that £ is a differential operator and construct its 
Green's distribution. The Green's distribution Q(x,t,y) is given by the semigroup acting on the Dirac delta 
function centered at y: 

g{x, t, y) = (e*^J(. - y)) (^) = ^ ^ e*^ ([A - - y)) (x) dA ^ ^ e*^G(x, y, A) dA. (2.1) 

Using variation of constants, solutions to the nonlinear equation can then be written as 

v{x,t)= / g{x,t,y)vo{y)dy ~ / / gy{x,t - s,y)Q{v,Vx){y, s) dy ds. (2.2) 
Jr Jo Jr 



To derive decay estimates for the Green's distribution, and hence for the integrals in (2.2 1, we can now deform 



the contour F in the rightmost integral in (2.1) pointwise for each {x,y). The key is that the resolvent kernel 



G{x,y,X) satisfies the ordinary differential equation 

AG - 6ix -y) = C - [fu{u{x))G]x, (2.3) 

so that ODE techniques such as the Gap Lemma [U [5] can be used to extend the resolvent kernel meromorphically 
across A = for each fixed (a;, y). It is then possible, see [71[33], to move the contour F for each (x, y) to extract 
the leading order behavior of the Green's distribution Q{x,t,y). 

We now illustrate the outcome of extending the resolvent kernel meromorphically across the origin A = 0. 
First, A = is a simple embedded eigenvalue with the exponentially localized eigenfunction u^. The associated 
adjoint eigenfunction is a constant V'l G I^^- Eigenvalues correspond to poles of the resolvent, and we therefore 
expect to obtain the term \^^Ux{x){4'i, •), which involves the spectral projection, when we extend the resolvent 
kernel across the origin. A second contribution should appear due to the essential spectrum: To determine its 
contribution, we consider the constant-coefficient equation 

\G-Gx. + A{x)Gx^5{x-y), A{x)^\ ^;^''~\ (2.4) 

Ju{u+) x>Q 
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for y < 0, say, where we neglected the shock profile, which we already accounted for through the spectral 
projection. To construct bounded solutions of this equation, we need to find the roots of the characteristic 
equation det(A — + fu{u±)v) = 0, which are given hy u k. and v « — A/a^ for A near zero. The latter roots 



are the dangerous ones as they give small exponential rates for the associated solutions of (2.4 1, and thus they 



will determine the next term, after the spectral projection, in the expansion of the resolvent kernel. In addition, 
we need to select those roots that give exponential decay of the resolvent kernel in (x, y) when A > 0. Using this 
information and the notation = {a^^^Ta~„ 



G{x,y,\) « -u^{x){i)i,-) + < 



, uq^j}, the expansion 

C^g--^2:/a+^j+Aiy/a" 

out. 
in ' 



for X < y <Q 
for y < X <Q 
for y <Q < X 



(1) 
(2) 
(3) 



for the resolvent kernel can be derived when y <Q for appropriate coefficients Cj that depend on the summation 
indices. The different cases in the equation above account for the transport of the perturbation along different 
characteristics, which we will discuss in more detail below. The inverse Laplace transform of the resolvent kernel 
is given by 



Q{x,t,y) 



+ E - 



E 



E 



C2 

-e 



/Ant 



C3 



Ant 



Thus, the spectral projection onto the cigcnfunction Ux appears, while the essential spectrum leads to Gaussians 
that move along the characteristics as they decay. It turns out that it is advantageous to use the term 



^(a;)E' 



in place of •), which uses an expansion of the adjoint eigenfunction in terms of the weak spatial eigenvalues 

V = —X/a^^. This leads to the decomposition G{x, t, y) — £{x, t, y) + G{x, t, y) of the Green's distribution with 



£{x,t,y) 
G{x,t,y) 



Ux{x)TT{y,t) 
j errfn 



E 



Cl 



y + «in^ 




E 



C3 

e 

/47rt 



t(t-l!//«" 1))^ 



''in. "out 



for y < 0, where we again use the notation a^^^ and a^^ to denote outgoing and incoming characteristics to aid 
in the following intuitive explanation of the above representation. 

The term £ consists of the eigenfunction and the inverse Laplace transform 7r(y, t) of the expansion of its 
adjoint eigenfunction. Hence, £ is an expansion of the spectral projection. The sum in 7r(y, t) is taken over 
the incoming directions, and it therefore tracks the initial data and nonlinear interactions resulting from the 
perturbation as they move in towards the shock and records their effect on the shock location. A sketch of the 
function tt is given in Figure |2] notice that Tr{y,t) 1 pointwise as t — *■ oo, and so UxT^ — > Ux{ipi, ■)■ The three 
pieces of G record, respectively, how the perturbation moves along the different characteristics that transport it 
in three distinct ways, as illustrated in Figure [2j outwards away from the shock (1); inwards towards the shock 
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Figure 2: On the left, a sketch of the expansion ^{y, t) of the adjoint eigenfunction is shown. On the right, the 
correspondence between the three terms in Q and transport along the characteristics is illustrated. 

and then reflected back out again (2); and in towards the shock, through it, and outwards on the other side of 
the shock (3). 

We can now outhne the nonhnear stabihty argument, which uses the above decomposition in the integral equation 



(2.2 1. More precisely, to track the shock location and remove the neutral direction Ux along the shock translates 



from (2.2 1, we write solutions as 

u{x + q{t),t) — u{x) + v{x, t). 
Exploiting the decomposition of the Green's distribution, we define the phase shift q(t) implicitly using 



'^{y,t)vo{y) dy- 



TTy{y,t- s) [Q{v,Vx) + qv] {y,s) dyds 



and find after some algebra that the perturbation v satisfies the integral equation 



v{x,t)^ / g{x,t;y)vo{y)dy 



Gyix, t- s;y) [Q{v, v^) + qv] {y, s) dyds. 



(2.5) 



(2.6) 



Note that the evolution of v is now governed only by the decaying part G of the Green's distribution. In [7l 133] . 
it was shown that pointwise bounds, which result from the above formulas for £ and Q, can now be used to 
establish existence of solutions to these integral equations and to prove that q(t) converges to a limit g*, while 
V decays algebraically to zero. 



2.2 Nonlinear stability of time-periodic shocks 

Given the success of pointwise estimates in establishing stability of stationary viscous shocks, we would like to 
use the same technique in the case of time-periodic shocks. For this approach to work, we must show that the 



resolvent kernel and the Green's distribution in equation (2.11 all have well-defined counterparts for a time- 
periodic linear operator. 

For time-periodic linear operators, the appropriate notion of the spectrum is the Floquet spectrum, and Floquet 
exponents a and the associated Floquet eigenfunctions u{x,t) are found as solutions to the linearized equation 

(ju + ut = Uxx-[fu{u{x,t))u]x, (2.7) 

where any eigenfunction must be localized in space and satisfy u{x,t + 27r) = u{x,t) for all t. Due to the 
nonuniqueness of Floquet exponents, it suffices to consider only a with — 5 < Imcr < ^ (see ^3.2 below). Based 



on (2.7 1 and comparing with (2.3 1, the resolvent kernel in the time-periodic setting should satisfy 



cru + ut- S{x - y)5{t - s) = u^x - [fu{u{x, t))u]x, 



(2.8) 



where the additional temporal delta function sets the initial time t^s. In order to solve (2.8 1, we first focus on 



(2.7 1 and write it as the first-order system 



^ -uf" 1 ,)-AU=:A{x,a)U 

\Ot + <7 + iuu{u)[Ux,-\ Ju{uj I 
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in the evolution variable x, posed on the space H™~^i{S^) x H"^{S^) for some to > 0. When this spatial 
dynamical system has an exponential dichotomy, given by solution operators <l>*'(a;, y, a) for x > y and (f>"(a;, y, cr) 
for X < y with <I>''(a;, x, a) + <l>"(a;, x, a) = id, then the inhomogeneous equation 



U,,^ A{x,a)U + H{x) 
can be solved uniquely, and the solution is given by 



U{x) 



<^%x,z,a)H{z)dz ~ / $"(x,z,cr)H(z)dz 



(2.9) 



(2.10) 



Thus, solving equation (2.8 1 is equivalent to using H{x) — (0, —d{x — y)S{- — s))'^ , which leads to 

-Pi$"(a;,y,a) 
Pi^^{x,y,a) 



G{x,y,a,s) = < 








X >y 
y > X 



as the resolvent kernel of ( 2.7 1 , where Pi projects onto the first component. Taking the inverse Laplace transform, 
we find the Green's distribution via 



g{x,t,y,s) = 



27ri 



'[G{x,y,a, smda. 



There are various issues that need to be addressed to make this approach work. Foremost among these issues is 
the regularity of the resolvent kernel G{x, y, a, s), since it is obtained by taking [H{x)]{t) = (0, ~S{x — y)S{t — s))'^, 
for which we cannot solve (2.9 1 in iJ™+3 x H™. 



3 Construction of the Green's distribution 

Consider the linearization 



Ut = u^x - [fu{u{x,t))u]^ (3.1) 



about the shock profile u{x,t). We say that Q{x,t;y,s) is the Green's distribution of (3.11 if, for each given 
Uf) e L-'^(M), the function 

u{x,t) ^ / g{x,t;y,s)uo{y)dy 



is a classical solution of (3.11 for t > s, and we have u{x, t) — > ^0(2;) as t \ s for almost every x € 



In this section, we shall show that (3.1 1 has a Green's distribution Q{x, t; y, s). However, knowing its mere exis- 
tence is not sufficient: to prove linear or nonlinear stability of time-periodic viscous shocks, we need to establish 
appropriate pointwise bounds for the Green's distribution. Thus, we shall construct the Green's distribution 
in a way that allows us to derive such bounds. Our strategy for finding the Green's distribution is as follows. 
Starting with the Green's function QQ(x,t) of the damped heat equation, which satisfies 

[dt-dl + l]go = 0, GoUs^Six^y), 

we will iteratively construct a sequence Gj of functions that satisfy 

[dt-dl + i]gj^[i~dx-fu iu)]g,-i , g, | = o. 

We shall see that will be become more regular as j increases and, in addition, the difference C/* of the sum 
g — X]j=o these functions and the desired Green's distribution g will become smoother as well. For a 
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sufficiently large we can then construct the difference as the inverse Laplace transform of a function G*, 



which, in turn, is given via the expression (2.101 for some regular Tl. Pointwise bounds for Q can now be derived 
immediately, upon exploiting their explicit construction. As we shall see in ^ we can also obtain pointwise 
bounds for by using spatial dynamics. The following theorem summarizes the existence result for Q and the 
pointwise bounds for Q. 



Theorem 4 Equation (3.1) has a Green's distribution Q{x,t;y, s), which lies inC'^OC} fort > s and is bounded 
uniformly in y for a; G M and t > s. The Green's distribution Q can be written as Q — Q + Q^, so that the following 
is true: For each a G M, there is a constant M > I so that the function Q obeys the pointwise estimate 

\d^g{x,t-y,s)\ < C(t-s)^^e-'"'""n'-T" , a;,2/eE, t>s (3.2) 

for Q < a <2. Furthermore, the function Q.^, is given as the contour integral 

g^{x,t;y,s) ^ — e''^[G^{x,y,a;s)]{t)da, for fixed > (3.3) 

2ni /„ i 

for a function G* that is defined pointwise in (x, t; y, s) and is analytic in a for a to the right of the Floquet 
spectrum S. 



In particular, the estimate (3.2) shows that, in Theoremjs] the term Q can be subsumed into the remainder term 



Q. Thus, once the preceding result is proved, it remains to establish pointwise bounds for Q^, to complete the 
proof of Theorem |3] this will be accomplished in Q-SjSj In the remainder of this section, we will prove Theorem |4] 
and establish additional properties of Qj and ^y*. 

3.1 Construction of 

Let Qo{x,t;y, s) be the Green's distribution 

Goix, t; y, s) = ] e-^"^*-^) (3.4) 
V47r(t - s) 

associated with the heat equation 

[dt -31 + l\Qo{x,t;y,s) = 0, t/o(a;, s; y, s) = 5{x - y), 
where x, y G M and t > s. We use Qq to define the functions Qj{x, t; y, s) recursively for j > I as solutions of 

[dt -dl + l]gj - [1 - a, • f^{u{x,t))]g,.u Gji^, s; y, s) = 0. 
Note that these functions are, by Duhamel's formula, given explicitly by 

gj{x,t]y,s) = / / t;o(a;,i;y, s)[l - % • /„(u(y,s))]^j_i(y,s;y,s)dyds (3.5) 



gQ{x,t\y,s)gj_i{y,s;y,s)dyds+ I I dygo{x,t; y, s)fu{u{y, s))gj^i{y, s; y, s) dy ds. 

We are interested in finding the Green's distribution g of ut ~ u^x ~ [/ji('")'"]k, which satisfies 

[dt -dl + dx- fu{u)]g = 0, g{x, s- y, s) = 5{x - y). 

Thus, seeking g in the form g = g^ + 5I]j=o that g is the desired Green's distribution if and only if 

g^{x,t;y, s) satisfies 

[dt-dl + dx- fu{u)]g,^[l-dx- fuiu)]gi, g,(x,s;y,s) - O. (3.6) 

We postpone the discussion of g^ to the next section and focus in the remainder of this section on pointwise 
bounds for the contributions C/j. 
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Lemma 3.1 Assume that f £ , then, for each j > 0, there is a constant C > 1 such that 



dlydtg,{x, t- y, s) < C{t ~ s)^-V^-l/5le-&-^ 
for multi-indices (a,/3) with < max{|Q;a;| + \(3t\, \oty\ + <k—\, and 

(9. + dy)\dt + dsfdlX^g, (x, t- y,s)\< C{t - s) ^ 
for < max{|a,| + \ay\ + + I7I + |<5| < fc - 1. 



I/3L- 



e C(t-a) 



(3.7) 



(3.8) 



Proof. Equation (3.4 1 can be used to obtain (3.7)-(3.8l for j = by direct computation. We can therefore 
proceed by induction to obtain estimates for j > I. The semigroup property for the damped heat equation 
Cut = Uxx — u imphes that 

, 1 x-hI^ t-I 1 \y-v\^ f-s ^ 1 [x-hI^ t-a 

(t-s)"5e c (g_s)-2e c(s-a) c dy = C(i - s)"5e ^(t-e) c , 



Using the induction hypothesis and (3.5 1, we can therefore estimate 
/■* 

\Gj{x,t;y,s)\ < 



{\Gq{x, t; y, s)||^j_i(y, s; y,s)\ + IdyGoix, t; y, s)||/„(u(y, s))^j_i(y, s; y, s)\^ dyds 



dS 



,-1 U^-i/l^ t-3 

< C{t-s)~e c 



which verifies the estimate (3.7 1 for a = /? = 0. 



To obtain (3.8) for a — f3 — 0, we observe that {dx + dy)Go{x, t; y, s) — {dt + ds)Go{x, t; y^s) — since Qq depends 
on its arguments through x — y and t — s only. Thus, after integration by parts in (y, s), we obtain 

ft 



{dx + dy)Qj{x,t]y,s) 



Qq{x, t; y, s){dy + dy)Q.j-i{y, s; y, s) dyds 



+ / c)ygo{x,t;y,s)fu{u{y,s)){dy + dy)gj^i{y,s;y,s)dyds 
Js Jm. 

+ / dygo{x,t;y,s)[dyfu{u{y,s))]gj-i{y,s;y,s)dyds, 

and similarly for (dt + ds). Using the induction hypothesis, we can now estimate these integrals as above to 
obtain ( |3.8[ ) for a = /3 = 0. Notice that we lose one degree of regularity for /„ S C''~^ for each application of 
(dx + dy) or {dt + ds), hence I7I + \6\ in total. 



Finally, using (3.8 1 to shift (y, s) to (y, s) derivatives where needed, and arguing by induction on j, \a\, and 
I7I, 1^1 , we may estimate 

m,ydl,g,ix,t;y,s)\< 

dt^df^goix, t; y, ~s) I |g,_i(y, ~s; y, s)| + \dl"^df^dygo{x, t; y, S)| |/„(w(y, S))e,_i(y, S; y, s)\) dyds 
{\go{x,t;yrs)\ \dl"^di^^gj-i{y,S;y,s)\ + \dygo{x,t;y,~s)\ |/„(S(y, S))^!"!^^ 5,_i(y, s; y, s)|) dyds 
plus terms that contain lower-order derivatives of gj, {dx + dy)gj, or (dg + dt)gj, yielding 



KydtGjMy,s)\ < Cit-syh^ 



C(t-s) c 



(t^rs) 



l + {t-s) 3 (s-s)^ds 



+C(i-s)-5e c / (S-s) 2 [l + {t-S)-^j 



dS 



< C(t-s) 5 e <=■(*-») c 



12 



as claimed, which verifies (3.7 1. A similar argument yields (3.8 1. Noting that each shift of x,t or y,s derivative 
costs one degree of regularity for /„ e C^~^ , and that we must shift \ay \ + la^l derivatives in the first integral 
and \oix \ + la* I in the second, we obtain the stated range of indices. ■ 



The computations in the proof of Lemma 3.1 may be recognized as parametrix estimates as in standard short-time 
parabolic theory [3l |5j [10] . 

The bounds for Q = X]^=o asserted in Theorem |4j are now a consequence of Lemma 
a e M and C > 0, there is a constant M > 1, chosen sufficiently large, so that 



3.1 



Indeed, for each fixed 



e C(t-s) C < 7\^e M(t-a) 



for all x, y e M and t > s. To complete the proof of Theorem |4] it therefore remains to verify the assertions 
about G*- 



3.2 Construction of Q^: 

In this section, we verify the assertions made in Theorem |4] about the contribution Q^, to the Green's distribution 



G and show that it is given by (3.3) for an appropriate function G«, which we shall refer to as the resolvent 
kernel. 



Recall that needs to satisfy equation (3.6 1, given by 



[dt -dl + dxfu{u{x,t))]g^{x,t; y, s) = [1 - dxfu{u{x,t))]gg{x,t;y,s), g^{x, s; y, s) = 0. 



(3.9) 



We shall use the Laplace transform to solve (3.9 1, since this approach will allow us to reformulate (3.9) as a 



spatial dynamical system, which facilitates the verification of pointwise bounds. 

It is convenient to shift the time-dependent coefficients, and we therefore consider the equivalent system 

[dt-dl + dxfu{u{x,t + s))]g4x,t)^ge{x,t;y,s), e*(x,0) = 0, (3.10) 
where we omit the arguments (y, s) in the notation for g^,{x, t) and use the notation 

ge{x,t;y,s) := [1 - dxfu{u{x,t+s))]ge{x,t+s; y, s). 
The term fu{u{x, t -\- s)) is smooth and 27r-periodic in t and can therefore be represented by its Fourier series 

/„(72(x,t + s))=^A.(x)e''=(*+^). 



Recall that the Laplace transform is defined by 

v{x, A) 



fcez 



v{x, t) dt 



for all A G C with ReA > /i, where n is chosen so that the integral is convergent. Taking the Laplace transform 



of equation (3.10), we obtain 



\g,{x,X)^dlg^,{x,\)^dx [^e'^^/,(a;)^,(a;,A-i^-) ) 



t{x,\;y,s). 



The effect of the time-periodic coefficients is that the equations for different values of A couple. Note, however, 
that the equations for CJ* at A = Ai and A = A2 couple only if Ai — A2 G iL. To exploit this fact, we define 



A = 



1 ^ 1 
-2<Ima<-, 
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which is motivated by |5] and similar to a Fourier-Bloch wave decomposition of periodic functions. For each 
(J, we can view the above equation as a system of infinitely many second-order ODEs in a;. In particular, if we 
define S^{x, a) := S*{x, a + in), and similarly for gi, we arrive at the system 



(a + in)g: = dlo: - [Yl ) 

VfeGZ / 



+ 5?(a;,tT;2/,s). 



(3.11) 



The following result, which we will prove in ^3.3 below, implies the existence of solutions to (3.111. Recall the 



definition (1.4) of the Floquet spectrum E of (3.101 



Proposition 3.2 Fix £ > S. For each a to the right of the Floquet spectrum S, the system 

aG, + dtG, = dlG, - [/„(w(x, t + s))G,], + ^ e'^'g^x, a; y, s), G^x, 0; y, s,a)^0 (3.12) 

riGZ 

has a unique solution G'*(x, t; y, s, a) in G^ that is 2tt -periodic in t. Furthermore, G^{x, t; y, s, a) is analytic 
in a and lies in Cy D Gl with respect to {y, s). 

Writing G^,{x, t; y, s, a) as the Fourier series 

G^{x,t;y,s,a) = ^ e'"*^"(a;, ct; s), 



we see that the Fourier coefficients Ql^{x,a;y, s) satisfy (3.11). Furthermore, the inverse Laplace transform 



1 /-Ai+ioo 1 /-M+s 

g*{x,t;y,s)^ — e^*g,(x,A)dA= — / e^*G.(x, t; y, s, a) da (3.13) 



of g*{x,X) is well defined for each fixed /i > and satisfies equation (3.10). This completes the proof of 



Theorem |4] subject to the proof of Proposition |3.2| which we will give in §3. 3| below. 



From independence of ( 3.13 ) with respect to /i > together with analyticity of G.^, with respect to a on {Re cr > 0}, 



we may conclude from Cauchy's integral theorem that 

e("+'/2)*a {x, t;y,s,a + i/2) - e^^-'^^^'G, {x, t; y,s,a~ i/2) 



(3.14) 



for all t. This is evident at a formal level, as is the more general property that e°'*G'* is periodic in a with period 
i or, equivalently, e~^'^^G^,{x,t;y, s, a) — G^,{x,t;y, s, a + ik) for all fc G Z, which follows since the left-hand side 



satisfies (3.12) with a replaced by cr — ifc. 



3.3 Construction of the resolvent kernel G^, 

In this section, we prove Proposition |3.2[ Throughout this section, cr will lie in the set f2, which denotes the 
connected component of C \ E that contains a — oo. We need to construct a solution G.^,{x, t, cr) of 

aG^+dtG^ = dlG^ - [fu{u{x,t + s))G^]^ + gi{x,t,a;y,s), G^{x,Q,(7) = 

that is analytic in ct e and 27r-periodic in t, where ge.{x, t, a; y, s) denotes the Fourier series of {g"(a;, cr; y, s)}. 
We rewrite this equation as the spatial dynamical system 

U, = A{x,a)U + ^,{x):=( ^ , ,!-,V+( ~, I' J (3-15) 

\Ot + a + }uu{u)[u^,-\ Ju{u) J y-gi{x,t,a;y,s) J 

in the evolution variable x, with u — u{x, t + s), where U{x) — {u, Ux)^ is a 27r-periodic function in time for each 
fixed X. 
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To construct solutions of (3.151, we first focus on the associated homogeneous system 



\at + (T + Juu{uj[Ux,-\ Ju{u) I 



(3.16) 



We consider the spatial dynamical system (3.161 on the Hilbert space K,„ = iJ™+2 (5^) x H™{S^), where m > 



will be chosen later. It can be shown that the operator on the right-hand side of (3.161 is closed and densely 
defined with domain y,„ i i; see [52]. Equation (3.161 is ill-posed, in the sense that solutions to arbitrary initial 



data in Ym may not exist: indeed, the leading-order operator 

(a. o) 

has spectrum given by {±\/i/c : A: G Z} and therefore does not generate a semigroup on Ym- Nevertheless, 



equation (3.16 ) provides a useful framework for analyzing the PDE (3.10 1, since it admits exponential dichotomies 



whose properties can be related to spectral properties of (3.101 



or M. Equation (3.16) is said to have an exponential dichotomy 



Definition 3.3 ([HI §2.1]) Let J = M+, 

on J if there exist positive constants K and k'^ < < and two strongly continuous families of hounded 
operators '^^{x, z) and <i>"(a;, z) on Ym, defined respectively for x > z and x < z, such that 

sup e-«^(--^)|i<i>^(x,z)|U(y„)+ sup e-''"(^-^)||<i>"(a:,z)|U(y„) <if, 
the operators P'^(x) := <i>^(x,a;) and P"(a;) :— ^^(x,x) are complementary projections for all x Cz J , and the 



functions ^^{x, z)Uo and $"(2;, z)Uo satisfy { 3. lay for x > z and x < z, respectively, with values in Ym for each 
fixed Ua G Y„^. 



It follows from [22l Remark 2.5 and Theorem 2.6] that (3.16) has an exponential dichotomy on M for each cr e il 



Furthermore, (TTJ Proof of Theorem 1] implies that the operators ^^{x, z, a) and ^^{x, 2, a) are analytic in a for 
fj G f2 as functions into L{Ym)- As a consequence of ^22,, §6.1], we can then solve the inhomogeneous equation 

= A{x,a)U + n{x) 

uniquely for each Ti G i'^(M, Ym) via the variation-of-constants formula 



U{x) 



and the solution satisfies 



<i>''{x,z,a)n{z)dz+ / ^''{x,z,a)H{z)dz 

-00 J 00 



U e H\R,Ym) n L^{R,Ym+i). 



(3.17) 



Our goal is to apply these results to equation (3.151, 



U^ = Aix,a)U + Ai{x), Ai{x):={ , 

\-9i(x,t,a;y,s}J 

for (T G ri. To establish the regularity of the right-hand side Ai{x), we repeat the iterative construction of the 



components Gj from ^3.1 for their Laplace-Fourier transforms, which is akin to the bootstrapping arguments 



carried out in §5.3] and Thus, we are led to consider the equation 

V,^Ao{a)V, Ao{a) 



r 

dt + (T + l 0, 
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which corresponds to the heat equation ut — u^x — u that we utihzed in ^3.1 Writing this equation in terms of 
its Fourier modes V", we obtain the system 







1 



which admits the exponential dichotomy 
$o'"(a;,?/,cr) 



in + a + l , 



yn 



1 



2\/a + 1 + in \ cr + 1 + in \/cr + 1 + in , 
1 / Vo" + 1 + m -1 



Vc+l+inla-j/l 



- V(T+T+in I a: — y I 



(3.18) 



2VaH^T+ln \ -((T+ 1 +in) \/(TTTTTn / 
In line with the definition of Qq as the Green's function of Ut = u^x — u, we consider the equation 



dxVo = Ao{(y)Vo + Ao{x), Ao{x) 





-5{x-y)6{t)^ 



(3.19) 



Using (3.17) and (3.181, the Fourier modes Vq of the solution Vo{x) are then given by 

$o"(a;,z, cr) 

\ --ya+l+inlx-yl 



1 



1 



+ 1 + in \ — sgn(a: — y)\/o' + 1 + in / 
The function Vi that corresponds to Gi can be found by solving the equation 

dxVi = Aoia)Vi + B{x)Vq, 

where 

B{x) := A{x, a) - Ao{a) = 





/««(w)[ui:, •] - 1 /«(""), 



(3.20) 



The equation for Vi can again be solved explicitly using the exponential dichotomy (3.18) for the Fourier modes 



To do so, let {ttk} and {bk} denote the Fourier components of fu{u) and fuuiu)[ux, •], respectively, (suppressing 
the dependence on s) and note that 

[m7o]»(x,y,a) = Ve-^^™l--^l f bn-k{x) V ,a„-.(x)|, 



where we set bo := bo — 1 and dropped the tilde. Applying the exponential dichotomy (3.18), we obtain 



(3.21) 



-^—y/o--\-l-\-i7l\x — z\ 



E 

fcez 



g-Vo"+i+ifc|z-y| 



v 



bn-k{z) sgn(a; ~ 7/)a„_fc(z) \ 



+ 1 + inV (T + 1 + ifc v^oH^X+ln 
= + sgn(a; - y)a„_fe(z) 



dz. 



Vo- + 1 + i/j 

We record the following estimate on Vi, which will be used below. By estimating the "worst" term in , we 
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obtain 



iV'il 



L2(R,L2(Si)2) 



J2 I |l>i"(x,y,a)pdx 



< 



< 



< 



< 



< 



n 
n 

^E 



k 

2 



^a„_,(.)e-^'™^l-«l 



L2( 

2 



|„-Vo-+l+ifcH II 



^E | a fEi 



|(7 + l + ifc|2 

We now reorder the second sum by defining j := n — k and use the fact that 

(i+inr)<c(i+|n-in(i+iin 



to obtain 



ll^^ill 



1,2(18,1,2(51)2) 



< 



CE r IE 



a,- z,o 



1 



\a + 1 + i{n - 



< ^lElKIU- (l^ + i|^ + b1^) ^ 

" + l + |a+l|5(|a+l|5 + |n|5) 

cE^^. 



< 



which is sufficient for convergence and guarantees that Vi G -L^(M, (L^(S'^))^). In the above derivation, we 
used the assumptions that / and u arc sufficiently smooth so that {afe} G £^'"(L!^(M)) with a = |, where 
i^''^{Lf{R)) = {{uk} : Y.k \k\"\\uk\\L^ < oo}. Note that Vi is analytic in a for Re<j > -1. 

From this point onwards, we can obtain the Fourier-Laplace transforms Vj of Qj inductively by solving 

In fact, the following lemma will allow us to repeat the bootstrapping procedure until we obtain any degree of 
smoothness in the inhomogeneity that we like, subject to restrictions only from the smoothness of f{u). 

Lemma 3.4 Assume that {ck}kez G ^^'"{L'^i^)) for some a> P > 0. IfV satisfies 

C 



[V{x,ym = Ee-*K(a;,y), ||K(-,y)|Ui < - 
uniformly in y, and V is defined by 



[V{x,y)\{t) = Ee-*K(a;,y), := / e-V-+i+in|x-z| ^ Cn-k{z)Vk{z, y) Az 
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then the functions Vn satisfy 



\\Vn{;y)\\h < 



C 



and 



C 



l + |n| 5+2/3 

where the constant C is independent of a and n. 

Proof. The assertion follows from a calculation similar to the one used above to bound V". We have 

2 



llK(-,y)| 



< C 



Vfeez / 



-V<T+l+inM 



L2 
2 



ke2 



LI 



< 



< 



< 



c 



- l + |„|H2/3- 

An analogous estimate holds for the norm of y). I 

The following lemma illustrates how the regularity of the functions Vj increases after each iteration of the 
bootstrapping procedure. 

Lemma 3.5 Pick any sufficiently small e > 0, then, for each j > 0, there is an rj > Q so that the function Vj is 
analytic for Recr > — \ and 



where :== H^{S^) for all k. 

Proof. The claims about Vq and Vi follows from their explicit formulas and estimates of the type given 
above: If we restrict cr to Recr > —i, we can also extract a factor e"'''^^^' from the convolution integral 



(3.211. Furthermore, applying dy and integrating by parts proves the claim about {dx + dy)Vi as we have 



||a,S(z)|U(y.„) <e-«l-l. 

Next, note that the coefficients {a^} and {hj} that denote the Fourier components of fu{u) and fuu{u)[ux, •], 
respectively, satisfy the hypothesis on {cj} in Lemma 3.4 for a value of a that is determined by the smoothness 
of the nonlinearity f{u). To estimate Vj, we use the fact that 



<i>l{x, z,a)B{z)Vj-i{z,y) dz + / ^"^{x, z,a)B{z)Vj^i{z,y) dz 



and denote by V" the least well-behaved of the two components of each Fourier mode V^". If V}_i satisfies 
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then Lemma 3.4 implies that 



illLi(R,L2(si)xL2(5i)) 



< 



R I 1 



^an-k{z)V'p{z,y) dz 



and 



j\\L^{M,L^(S^)xL^(S^)) 



< 



< 



< 



^E 

ra6Z 

^E 



l + |n|J+3 



I 1 



fcez 



L2( 



which gives the LiJ. and bounds on Vj. Note that we can again extract a factor e^^'l^^^' from the convolution 
integrals. To address the iJ^ bound, note that the derivative dx falls only on the exponential Q-V<7+T+m\x-z\^ 
and the resulting extra factor of \n\i leads to the loss of smoothness in time stated in the lemma. ■ 



We now turn to equation (3.15), which is equivalent to 

dxU = AU + B{x)Vi{x). 

If we choose £ = 3, then 

P2BV3 e L2(M,ij2+7(5'i)) ni/i(M,i/i+''(S'i)) 



(3.22) 



for each fixed < 7 ^ 1, where P2 denotes the projection onto the second component. Thus, if we set m = 2 + 7 
in the definition of the underlying space so that = H^+^{S'^) x H^+'f{S^), then BV3 € L'^{R,Y^). As 
discussed above, (3.221 then has a unique solution C/» — (G^,dxGi,), which depends analytically on a and lies in 
so that 

/ 



, dxG^ 



In particular, [G* {x, a)] (t) and [dxG^, (x, a)] (t) are continuous functions. Inspecting ( 3.22 ) and using the regularity 
of we find that 9^G* e H];.{R, H^'^'^{S^)) so that (3^G* is also continuous in {x,t). 

It remains to discuss regularity with respect to (y, s). Differentiability with respect to y is a consequence of the 
iteration scheme together with Lcmma |3.5[ Differentiability with respect to s follows similarly: since we replaced 
t by t + s at the beginning of our analysis, the initial inhomogeneity Ao(x) from ( |3.19 1 does not depend on s, and 
the dependence of Vq on s is only through B{x), that is, through fu{u{x, t + s)) and its derivatives. In particular, 
the recursive construction of Vj shows that dsVj lies in the same space as Vj. Once we shift back to the original 
time variable, s-derivatives of Vj become equivalent to i-derivatives and therefore lose one degree of regularity 
in time. 

This completes the proof of Proposition |3.2[ We record as a corollary the representation of G* given through 



the variation-of-constants formula (3.171: 

Corollary 3.6 The resolvent kernel G* can he represented as 
[G.{x,y,a-s)]{t) - Pi 



z, a)B{z)V3{z, y, cr; s) dz 
z, a)B{z)V3{z, y, a; s) dz 



it) 
it), 



(3.23) 



19 



where ^^'^ is the exponential dichotomy associated with the operator A on the space Ym, and Pi is the projection 
onto the first component. 

In ^ we will show that the exponential dichotomies <i>''(x, y, a) and <i>"(a::, y, a) can be extended meromorphically 
by separating out the translational and essential eigenmodes. The most natural way to transfer this result to 
the full resolvent kernel G{x, y, a, s) is via the formal expression 



G{x,y,a;s) = 



Pi<i>''{x,y,a)5{-) x>y 
-Pi$"(x,y,cr)(5(-) x<y, 



(3.24) 



where • denotes the argument t. However, as discussed above, we cannot apply the dichotomy directly to S{t). 
If one could argue, possibly using test functions and the uniqueness of strong solutions to equation ( 3.15[ ), that 
the Green's distributions given through G = G* + ^jGj ^^^d via (3.24) must be equivalent in a distributional 



sense, then one could instead work directly with the dichotomies via (3.24). However, we do not know how to 



make such an argument work because test functions would be smooth in t. In order to solve the initial value 



problem associated with (3.1), we must work with initial data of the form S(t)uo(x), which are not smooth in 



time. The function S(t) represents the spreading of the initial data amongst all Fourier modes and is a key aspect 
of the dynamics. This appears to be an important distinction between the time-periodic and time-independent 
problems. 



4 Meromorphic extension and bounds for the resolvent kernel 

The goal of this section is to extend the resolvent kernel G^{x,y,a; s) that we constructed in the last section 
meromorphically across the essential spectrum near a = and to derive sharp pointwise bounds for G* with re- 
spect to {x, y). To state our result, we define the spatial eigenvalues i^ficr) as the N solutions of the characteristic 
equation 

det(i^^ - fu{u±)v - cr) = 
that are close to zero when a is close to zero. These eigenvalues are analytic in a and have the expansion 

(a) = + + 0(|aH, J = 1, . . . , ^. 

We shall also use the notation and to denote the spatial eigenvalues associated with the outgoing and 
incoming characteristics a^^j and aj^, respectively; see Definition |l.l[ The following theorem is the main result 
of this section. 

Theorem 5 Assume that Hypothesis \ (H1 j] is met and that the shock profile u{x,t) is spectrally stable so that 
(S1)~(S4) in Definition 1.2 are met, then there exist positive constants C, rj and e so that the following is true: 
The resolvent kernel [Gi,{x,y,a; s)](t) has a meromorphic extension in a into {cr G C : Rea > — e} and can be 
written as 

G*{x,y,a; s) = Ei{x,y,a; s) + E2{x,y,a; s) + R{x,y,a; s), 
where the terms Ej have a pole at a — 0, while R is analytic in cr. For y < 0, we have 

E,ix,y,a,t;s) - -J2Mx,tKiniy.sfe"'^^^"^y (4.1) 



E2{x,y,a,t;s) = - V Mt(a;, i)?2 i„(y, s)^e 



for appropriate functions Ij^miu, s) that are 2'k -periodic in s, and a symmetric representation holds for y >0. The 
remainder term R satisfies the following pointwise bounds, where a is any multi-index with < la^;! -|- |q!j^| < 2 
and < \at\, \as\ < 1.' 
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(i) For X > y > 0, we have 



sup \R{x,y, a, t;s)\ < C ^ e'-'^'^'e"" 



S.t 

(a) For X > > y, we have 

sup \Rix,y,a,t;s)\ < C e^o-^-^^e-'^rJ-b, 

S,t _ 4- 



(Hi) For > X > y, we have 

sup \R{x,y, a, t;s)\ < C V e'^^f'^'^e-^inW^, 

s,i _ 



Symmetric bounds hold for x < y. 

In the remainder of this section, we prove Theorem [5] Due to Corollary |3.6[ it suffices to extend the two integral 
terms 

$'(x,z,cr)S(z)T4(z,2/,cr;s)dz+ / ^""{x, z,a)B{z)Ve{z,y,a; s) dz, (4.2) 



with ^ = 3, in the representation (3.231 of the resolvent kernel G*. Since B does not depend on a, and Ve 
is analytic near cr = 0, we need to extend the exponential dichotomies ^^{x,z,a) and ^^(x,z,a) for x ^ z 
with x, z G M. First, we use spatial dynamics to extend the exponential dichotomies on the half lines M"*" 
and E~ analytically across a — 0. Afterwards, we use the assumptions on the Floquet spectrum to construct 
a meromorphic extension of the exponential dichotomy on M and derive pointwise bounds for this extension. 



Finally, we transfer these bounds to the resolvent kernel G* by estimating the integrals (4.2) 



4.1 Analytic extension of the exponential dichotomies on 



Consider the spatial-dynamical system (3.161 



yot + (T + fuu{uj[urc,-\ iu{u)j 

on the space Y„i = x iJ™ for m > 2 fixed. For Rccr > 0, this equation possesses the exponential 

dichotomies <i>*'(x, cr) and <i>"(x, y, cr), which are defined and analytic in a for x > y and x < y, respectively, 
with x,y G M.. Our goal is to construct analytic extensions of these dichotomies separately for x,y G and 
x,y G K~ from Re cr > to a small ball B^{0) centered at cr = 0. Throughout this section, e denotes a positive, 
and possibly small, constant that we may adjust during the arguments to follow. 

First, we consider the asymptotic equations 

Ux=(^^^ J \u=:A±{a)U. (4.4) 

\Ot + (T Ju{u±) 



21 



-2 

X 

X 


X 


2 


• 


X 

X 




• 


• 

• 



X 


V 


2 

X 


• 


X 

X 


• 




• 

• 



RecT > 



Rccr < 



Figure 3: The spatial spectrum of A+{o-) is shown for Ke a > [left] and Retr < [right]. The spatial eigenvalues 
that are stable for Re a > are indicated by crosses, while unstable eigenvalues are shown as bullets. As indicated, 
the N small eigenvalues uj' cross through the imaginary axis as Re a changes sign. 



Using the fact that the operators A±{a) leave the 2iV-dimensional subspaces spanje'^^'F : V G C^^} C 
invariant for each fc S Z, it was shown in |24j that their spectrum is discrete and given by the spatial eigenvalues 



°' + ^-^VK^]' + 4('T + ifc), j-l,...,iV, fceZ, (4.5) 

where the are the nonzero, real, distinct eigenvalues of fu{u±) guaranteed by Hypothesis (HI) Furthermore, 
there is an 77 > so that these eigenvalues have distance 3?7 from the imaginary axis, uniformly in a G i?e(0), 
except for the N spatial eigenvalues 

-f(-) = -^+j|j^+0(|an, J- = l,...,iV, 



which arise from (4.5) when setting k — and expanding in a near zero. The eigenvectors of A±{a) associated 



with the eigenvalues vf{<j) do not depend on t and are given by 



1 




V±(a) I J rf , j = l,...,7V, 



where are the right eigenvectors of fu{u±) belonging to aj . Key to our analysis is the fact that these 
eigenvectors are linearly independent and analytic in a for all a near zero. 



For Re cr > 0, the eigenvalues v^{<^) move off the imaginary axis, and the operators A±{(j) are hyperbolic; see 
Figure [3] In fact, since a",^ < < ar and < < a^y^, we see that 

Rei'(^t(cr) > and Rei/^(.(cr) < for Recr > 0, 

so that these spatial eigenvalues contribute respectively to the unstable eigenspace £'"(cr) of A-{(t) and the 
stable eigenspace i?+(cr) of A+{(t) when Rcct > 0. Thus, we define the subspaces 

7^it(a) = span{V±(a) : af ^ 0}, 7^±(a) = span{V±(a) : af ^ 0} 

of outgoing and incoming modes, respectively, which are analytic in ct g B^(0). Similarly, we can define the 
spectral subspaces E±{a) and E}^{a) belonging to A±{a) that correspond to the stable and unstable eigenvalues 
with real part less than —877 and larger than 3r], respectively. Using these definitions, it follows from the above 
discussion and the results in [24 that the decompositions 

E^'^ia) ® 7^-,t(a) ®ET{a) ® 7^,^ (a) = Y^, E^_^{a) ® 7^+ ^(a) 0^7(^1) ® TZ+{a) = 

^ V ' ^ ' 

= E^{(y) for Recr>0 =E\(<j) for Recr>0 
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and the associated projections P"(cr) and P^{(j) exist and are analytic for a G i?e(0). In particular, the unstable 
subspace £'"(<j) of the stable subspace E^{<j) of A+, and their spectral complements can be extended 
analytically from Re cr > to the ball B^{0). 



Lemma 4.1 There are positive constants C and e so that {4--3) has an exponential dichotomy <I>^"(a;, y, cr) on 
that is analytic in a Cz -Bc(O) and satisfies 

||<i>;(x,y,a)|U(y„) < C^e''-('^)(--^), x>y>0 (4.6) 

\\<i>l{x,y,a)hiY^) < C5]e^i(^'(--^), y>x>0 

and 

||5,$;(x,y,a)|U(^,„,^^_^) < c(|a|+e-''l--^l)^e''o+u.(-)(-v)^ x>y>Q 

\\d^:,^l{x,y,a)\\Lf^Y^^Y_^) < C(|a|+e-''l--^l)^e-i(^)(--''), y > a; > 

for a G i?e(0) . Furthermore, the associated projection P^{x, a) := ^^_^{x, x, cr) on M+ satisfies P^{x, cr) P^{o') 
as a; — !■ cx). T/ie same statement with symmetric bounds holds for dichotomies on . 

This result can be viewed as an infinite-dimensional version of the Gap Lemma, which was established in finite 
dimensions in [H [5] . 

Proof. Since the distances of the strong stable and strong unstable spectrum of A+{(t) to the imaginary axis 
are larger than Srj, we know that the shock profile u{x,t) approaches its limit exponentially with rate Srj in 
the C^-norm as a; — > oo. 

We shall use the following result: If there are constants < such that the spectrum of the asymptotic 
operator A+{(j) is the union of two spectral sets defined by eigenvalues with real part respectively less than 



and larger than k", uniformly in cr S -Be(O), then (4.3 1 has exponential dichotomies on K+ with rates and k" 
as outlined in Definition |3.3| and these dichotomies can be chosen so that they are analytic in a. Furthermore, 
as X — > oo, the associated cc-dependent projections converge exponentially with rate min{3ry, |k" — k^\} to the 
spectral projections of -4+(cr) associated with the two spectral sets. This claim follows from [17, Theorem 1] 
upon using exponential weights, and we refer to |25j for further details. 

Choosing = —2r] < ~-rj = k", we find analytic dichotomies <i>^(a;, y, cr) and <I>™(x,y, cr) corresponding to 
solutions that decay with rate at least —2ri as x increases and solutions that grow not faster than with rate 
rj in backward time. Similarly, picking nf = r] < 2ri = k", we obtain dichotomies <i>^(a;, y, cr) and <i>""(a;, y, cr) 
that correspond to solutions which grow with rate at most as a; increases and solutions which decay with 
rate at least —2ri in backward time. The difference between these dichotomies is whether we subsume the N 
center directions that belong to the small spatial eigenvalues V^icr) into the unstable or the stable part of the 
spectrum. Using the convergence of the associated projections P'^{x,a) and P^{x,a), we see that the subspace 
E'\_{y,a) := RgP!j^{y,a) n Rg P™(y,CT) has dimension N for all y > and all cr. Furthermore, solutions U{x) 
with initial data U{y) in E^{y,a) exist for all x > with U{x) E E'^{x,(t), since we can use <i>^jf (x, y, cr) to 
evolve for x > y and <f>™(a;, y, cr) to evolve for x < y. Thus, we successfully isolated the A^-dimensional center 
directions from their infinite-dimensional stable and unstable counterparts on which we already have exponential 
dichotomies <i>'!f'(a;, y, cr) and <I>""(a;, y, cr) that are analytic in cr € -Be(O). 

Next, we decompose the A^-dimensional center space E'i^(y,a) into two complementary subspaces which are 
composed of solutions that converge with uniform rate 2ri to 7^out(''') oi' to ^hi(^)) respectively, as a; ^ oo. We 
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proceed as in §4.3]. First, we write (4.3) as 



^ [A+{<7) + B+{x)]U, B+{x)^ 











(4.7) 



where ||S+(a;)||i(y^) < Cc ■^''l^' for x > 0. Pick an index j so that > is an outgoing characteristic. We 



seek a solution U{x) of (4.7 1 of the form 

U{x) = e''/(")("-^)V+(a) + Vix), 



(4.8) 



where L > is some large constant and we require that |V^(-t)|y„ < Ce ^''^ as x oo. To construct V{x), we 
consider the integral equation 



Vix) 



^ J c-^+(")^+('^)("-")F;(a)S+(z) [F(z) + e''/(")(""^)V+(a) 



dz 
dz 



(4.9) 



for X > L. Upon fixing a sufficiently large L, it was shown in [20 , §4.3 and (4.12)] that (4.9 1 has a unique solution 
Vj{x, a) for X > L that grows with rate at most as a; — > oo. Once this is established, one can show that this 
solution depends analytically on a and, in fact, converges exponentially with rate 2r] to zero as x — > oo, because 
the first integral term becomes zero in this limit. In order to construct this solution for all a: > 0, we now need to 
flow it backward from a; = L to a; = 0. However, we cannot necessarily do so in this infinite dimensional setting 
because we do not know whether the initial data Vj{L, a) lies in E'^{L, a). However, it is easy to see that 

i?S,t(i,a) [RgP;^(i,a) ®span{y,(L,(7) : a+ > 0}]n El{L,a) 

has the same dimension as TZ^^^{a), possibly after making L larger. 

Proceeding in the same fashion for initial data in TZ^{a), we can construct an analytic complement Ef^^{L,(T) 
of £^Q^^(L, cr) in E'^{L,a). Since we can evolve initial data in E'^{L,a) for all a: > 0, we can define an analytic 
and invariant decomposition in E^{x,a) for each a; > 0. Adding <i>^(x, y, a) and <i>""(a:, y, cr) to the center 
evolutions that we just constructed defines an analytic extension of the exponential dichotomy on ^^_^(x,y, a) 
and (f>" (a;, ?/, cr) into Bg{0). The bounds stated in Lemma 
exponential bounds for Vj{x,a). 



4.1 



are a consequence of the ansatz (4.8 1 and the 



4.2 Meromorphic extension of the exponential dichotomy on 



We define -E+((t) and i?"(cr) to be the ranges of the projections P^(0,cr) := (0, 0, ct) and P"(0,cr) := 
<i>" (0, 0, cr), respectively. It follows from [17, Theorem 2] that the exponential dichotomies ^^_^{x,y,a) and 

separately on M+ and , fit together at x = y = to produce an 
E'^{a) — Yfn- Thus, if this equation were true for all a near 
for all such a, which could be constructed 



4.4 



^^{x,y, a), which we defined in Lemma 
exponential dichotomy on K if and only if E^ (ct) 



zero, then (4.3 1 would admit an analytic exponential dichotomy on 



explicitly through the analytic projection onto E^{a) with null space E^{a); see [17} (3.20)] and (4.24) below 



As we shall see in Lemma [4~2] below , the direct sum decomposition of Y„i through stable and unstable subspaces 
fails at (7 = 0, due to the presence of the embedded spatial and temporal translation eigenmodes. Therefore, we 
instead show that the projection onto E^{a) with null space i?"(cr) has a meromorphic extension in a, with a 
pole at (T = 0, which we can use to construct a meromorphic exponential dichotomy on JR. 



Consider (4.3) 



U,^Aix,a)U 





dt + a + fuu{u)[u^, ■ 




(4.10) 
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and its formally transposed equation 



-dt + a + fUu)[u,,-] 



1 



inn) 



W, 



taken with respect to the real inner product in X = L'^{S^) x ^^(S'l). A calculation [221 El] shows that 

d 



da; 



{W{x),U{x))x =Q Vx 



(4.11) 



(4.12) 



for solutions U{x) of (4.101 and W{x) of (4.11) 



We now set a = and consider the resulting equations 



\Ot + Juu{u)[Ux,-\ Ju{Uj I 



and 



W,^-A{x,OfW = - 



^0 -dt + fLmu.,-] 



1 



/J(") 



(4.13) 



(4.14) 



Equation (4.131 admits the two linearly independent solutions 



Ui{x) = d^{ , U2{x)=dti , 

which are defined for x e M and decay exponentially to zero as a; — > ±oo. Next, consider the adjoint equation 



(4.141, which admits the solutions 



W{x) 



(4.15) 



where w G M.^ is arbitrary. These solutions can be used to define N bounded and linearly independent solutions 
by substituting the N basis vectors ej in for w. As shown in 24J, these solutions originate from the smooth 
functional 



which is conserved under the evolution of the system 



27T 



[v-f{u)]dt, 



Vx 



dtu + fuiu)v. 



on K„i that time-periodic shock profiles with period 2n satisfy; see (1.2) 



As outlined in sjl] the condition (1.5) in Hypothesis (S3) implies that there is a unique nonzero vector ipi 



pN 



up to scalar multiples, that is perpendicular to the outgoing eigenvectors of fu{u±) so that ijii -L [Rtui ® ^outl- 
We define 



^1 



to be the associated solution of (4.14). Hypothesis (S4) implies that a second solution of the adjoint equation is 
given by 

^2 r 



^2{x) 



where ip2{x,t) appears in (S4). Assumption (1.6) implies that ^'i and ^'2 are linearly independent. Recall the 
definition i?±(o-) := RgP^(0, a) for j — s,u. 
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Lemma 4.2 The linear mapping 



is Fredholm with index zero for all a G -Be(O). Furthermore, we have 

EXiQ) n EliQ) = span{[7i(0),!72(0)}, [^^(0) + EliQ)]^ = span{^'i(0), ^-2(0)}, 
and E%{a) E'iia) = for all a G -Be(O) \ {0}. 

Proof. Spectral stability of the shock profile together with [5TJ §4] implies that is invertible, and therefore 
Fredholm with index zero, for Reu > 0. Furthermore, [TTJ Corollary 1] implies that the nuUspace of i(0) is 
finite-dimensional, while [T71 Comment on p. 273] and [22 Lemma 6.1 and §6.2] show that the range of t(0) is 
closed and has finite codimension. Thus, i{a) is Fredholm with index zero for cr = and hence for all a G i?e(0), 
possibly after making e smaller, as the set of Fredholm operators of a given index is open. 

Next, it is clear that Ui{Q) and [/2(0) lie in E^^{Q) n i?" (0), and Hypothesis (S2) implies that this space does not 



contain any other initial data that lead to nontrivial localized solutions of (4.131. Hence, the proof of Lemma 4.1 



implies that any other element in this intersection corresponds to a solution V{x) of (4.13) for which V{x) 



converges to a nonzero element of H-^^^ as a; ^ 00 or of H-^^^^ as x — s- — cx). Suppose the former case occurs with 



V{x) -> y+ = (r+,0)^ G 7^+ t \ {0} as x 



00. Since U^^^ n UZ^^ 



{0} by (S3) and because Ti-ont is invariant 



under fu{u±), we can pick a ■0+ € L-^^ so that {tp^, fu{u+)r+) = 1 and {ip^, fu{u^)r^) = for all r_ G 



Define the associated solution "^(x) of (4.14) through (4.15) and observe that {^>{x),V{x))x = —1 for all x by 
construction and (4.12). However, V{x) converges to R~^^^ (or to zero) as a; — * —00, and we reach a contradiction 
to our choice of ■)/'+• This proves our claim about £'+(0) n E^{0). A similar argument shows that ^i(O) and 
^"2(0) are perpendicular to the range of t(0) and therefore span the complement of the range as claimed. 

We can regard t{a) as being analytic in a by viewing the equivalent operator 



L{a) : EliO) X E'^iO) 



Y 



(1^-^ v^) ^~> p;(o, cr)y" - p'lio, (T)y". 



Thus, t{a) has a nontrivial nuUspace either for all a or else only for a discrete set of a in its region of analyticity. 
For cr > 0, a nontrivial nuUspace of L{a) corresponds to a Floquet eigenvalue, and Hypothesis (SI) precludes 
their existence. Hence, we conclude that (.(cr) is invertible for all a G B^{0) \ {0}, possibly after making e smaller. 



Lemma 4.2 implies that there are closed subspaces Eq, Eq, Eq^ , and Eq of Ym with 

E^„®E^®E^®E^ ^Yra, 

where 

E^' = EliO) n E1{0) = span{f7i(0), 1/2(0)}, 



(4.16) 



E^ = [El{0)+Eim^ = span{vI/i(0),vI/2(0)}. 



Note that E^ is not uniquely determined, and we shall use this freedom below in Lemma 4.3 to make a specific 



choice that simplifies the estimates. We define V to be the projection onto Eq Eq with null space Eq Eq 
Lemma 4.3 For each a G -Be(O) \ {0}, there is a unique mapping h^{a) : -E"((t) E^{a) so that E^{cr) 



pt 



graph /i+ (cr) . Furthermore, we can choose Eq subject to (4. 16) so that /i+(cr) can be written as (1 — V)h^(a)V , 
where 



h+{a):E^®El 



El®El\ 



/i+(a) = /i+(a) + /^+(a). 



26 



and h'^{<7) is analytic for a € B^{0), while /ip (cr) is given by 



h+{a){V\V^')^-{OjloV^), 
a 



where Mq : Eq — > Eq^ is invertible and has the matrix representation 



(■01,M)RiV 

(V'2,U2;)L2(5.i)da; / (7/;2,Mt)L2(5i)dx 



Mo = 

with respect to the bases {C/i(0), C/2(0)} and {^'i(O), ^'2(0)}. Similarly, for each a £ 3^(0) \ {0}, there is a 
unique mapping h^(a) : E^_{<j) E^{a) so that E^(a) — graph (cr) . This mapping has a meromorphic 
representation analogous to the one given above for h^{ij), but now involving the matrix ~Mq. 

In particular, ft.+ (cr) is meromorphic on i?e(0) with a simple pole at cr = 0. 

Proof. Our proof mimics Lyapunov-Schmidt reduction. We use the coordinates 

(ys^ VP\ y", v^) €E'o® E^' ®E^®E^ 

and indicate the range of operators by the appropriate superscript: the mapping g^'^{cr), for instance, maps into 
E^®Et. 

Lemma 4.2 implies that we can write i?" (cr) uniquely as a graph over 

£^0 ffi ^0* with values in El®E^. Thus, 

there are unique analytic mappings h^j{o') and h^{<^) with 

E'^ia) : V^V'' + VP' + /i"i(cr)y" + hl{a)VP' + hf{a)V'' + ht{cr)VP'. 
Let CT = 0. Setting = 0, we find that 

yp* + hi{o)vp' + ht{o)vp' e £;:^(o) V e E^\ 

Since E^^ C £^'^(0), we conclude that h^{0) ^ /i|(0) 0. Next, set V^p* = 0, then 

+ /i"i(o)i/" + /if (o)y" e E'^io) V g e'^. 

The only requirement for Eq is that £'q ® Eq — E^{0). Upon replacing Eq by graph /if (0) C i?"(0), we can 
therefore assume that /if (0) = 0. Thus, the above discussion shows that we can write E]!_(a) as 

E'Lia) : F = + yP* + a/i'^(CT)(y" + V^p*) + /i''(cr)y" + ah%a)VP\ 

where all mappings in the above expression are analytic in a. 
Since the subspaces E^{a) and E^{a) are analytic, we also have 

Elia): V = V' + VP' + ag'''^'{a){V'' + VP') (4.17) 

Elia): V^ = V^" + F'^ + cr/'P*((T)(y" + l/'^) 

for unique mappings g^"^ and g^'^^ that are analytic in a. 

We need to write E"i{a) as a graph over i?"(cr) with values in E^{a'). Thus, consider 

_^ ^pt _^ ^/jV'(cr)(t^" + v^pt) + h%a)V'' + ah%a)VP' (4.18) 
= [y" + + a/'P*(CT)(y" + V^)] + [V" + VP' + ag''^{(T){V' + V^')] , 
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where we need to express {V^, V^*^) in terms of {V^, V^) so that (4.18 1 is true. Upon writing (4.18 1 in components, 
we see that 



Substituting these expressions into the stable component of (4.18), we obtain 

h%a)[V'' + (Tg"((T)(y" + 1/P*)] + ah%a)[V'' + + ag''^P\a){V'' + V'^')] = + ag%a){V'' + V'''), 
which wc can solve for so that 

+ah''{aW + V^ + CT/'P'(a)(l^" + V^)] - <Tg'{a){V'' + V^)'^ 
=: hl{a)V'' + ahl{a){VP' + V^), 

where hj{(T) is analytic in a. Finally, the £'^-component of equation (4.181 is given by 

tT/i'''(cr)(y" + crg"((7)(/i^i(f7)F" + f7^((j)(yP' + V^) + yP') + V^' + a<?P*(a)(V^" + 
= + ag'^{a){hl{a)V'' + ahl{a){VP' + + VP'), 



(4.19) 



which is of the form 



id^^ + ahf{a) 



h^{a) - /(a) + ahtia)] V^' + <jht{<j)V\ 



where all mappings are analytic in a. Assume for the moment that 



M(0) = /i'^(O) ~5'^(0) : £;P* 



is invertible. We then have 



a 



where M(0) = [/i'^(O) - g'^iO)]-^, and M{a) and hl\a) are analytic in a. Using also ( |4.19| , we see that 

(V^, 1/Pt) = /i+(a)(F", V^) + -(0, M{0)V^) 

a 

as claimed. 

It remains to prove that M(0) = h'^{0) — g^'{0) is invertible and to derive an expression for its matrix representa- 
tion. Hence, we need to find expressions for the i?^-components of the graphs of E^{a) and E^{a) over Eq'. We 
start with El{a) = Rg $"^(0, 0, cr). Using the definition ( |4.17| of .g'''(cr) and recalling that X = L'^{S^) x L'^{S^), 
we see that 



(^',(0),/(0)[/,(0)>x = (*^(0),D„$;(0,0,a)U=oC/,(0)) 



X, 



(4.20) 



since the i?^-components of the projection P^{0,a) and the parametrization (4.171 differ only at order 0(cr'^). 
To derive an expression for Do-<i>^(0, 0, 0)C/j(0), we recall that ^^_^{x,0,(j)Uj{0) satisfies equation (4.3 1, which 
can be written as 

'0 0' 



A{x,0) + a 



1 0, 



U. 



Thus, V{x) = D^<i>^(a;,0,0)C/j (0) is a bounded solution of 





1 0, 



C/,(x), 
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where we used the fact that ^^_^_{x,0,0)Uj{0) — Uj{x). Using the construction of ^^_^_{x , y , a) in Lemma 
know that there are analytic coefficients ajk{cr) so that 



4.1 



we 



p-i 

<i>l{x, 0, a)C7, (0) = 0, a)[7, (0) + ^ aajk{cT)Vk{x, a), (4.21) 

k=l 

where the factor a arises since Uj{x) = <i>^(a;, 0, 0)C/j(0) decays exponentially with rate at least 877. Thus, V{x) 
is of the form 

V{x)^<P^4x,0,0)V,+ + y'%-(:.,z,0)r ^]u,{z)dz + J^'^^ix, z,0) T |Jjc7,(z)dz 



particular solution converges to as x — *oo 



for an appropriate Vj/": indeed, (4.21 ) shows that V{x) — + Ti-tut as x ^ oo, as does the particular solution in the 
above expression; thus, the only other contribution to V{x) can come from <i>^(a;, 0, 0). Hence, 



+ 

oo 



and, setting ui := and U2 := Ut, we obtain 



oo 


oo 



vl/,(0),<i>™(0,z,0) r ljU.iz)) dz 



ci>™(0,z,0)^vl/^(0), ( ^ JUj{z)) dz 



/ X 

'0 0' 
1 0, 



X 



I \U,{z)) dz 





] Uj[z) j 
(?/;,(z),Mj(z))i2(5i)dz, 



where we have used the fact that <i>™(0, z, 0)-'"^'i(0) = ^'i(z); see [251 Lemma 5.1]. Proceeding in the same 
fashion for E^{a), we find that 

_ .0 



and therefore 



M(0) = /l'^(O) -/(O) = ( r (V'.(z),U,(z))L2(5r)dz) . 



Furthermore, since ipi{x,t) — -0i £ K", we have 



27r /• / /•27r 

(V'i(a;,t),wt(a;,i))Riv didx = / (ipi, ut{x,t)dt) dx = 0, 



R ^0 



which proves that il/(0) is lower triangular. Hypotheses (S3) and (S4) imply that the diagonal entries of M(0) 
are nonzero, so that M(0) is invertible. 

This completes the proof for /i_|_((t). The proof for h^{a) is analogous, but we need to be careful with the signs: 
The integral representations of the solutions V{x) given above stay the same, but the roles of and g"^ are 
reversed. Thus, the matrix that appears in the representation of h~(a) is —Mq, and not Mq. ■ 
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We define 



<^>l{x,y,<j)Pl{y,<j) x>y>0 
{l-Pl{x,a))<i>l{x,y,a) y>x>0 



(4.22) 



and similarly for a;, y < 0. As in [T71 (3.20)], equation (4.221 defines an exponential dichotomy on M+ with 
projection P^{x,a), since RgP^{x,cr) — RgP^(a;, cr) and (1 — P^(x, cr))(l — P^{x,a)) = 1 — P^(a;,cr) for all 
a. Similarly, P"(a;, cr) is the projection for the exponential dichotomy ^t{x,y,a) and <E>l(x, y, cr) on M^. By 
construction 

Pl{0,a) = 1- P'^{0,a) Vctt^O, (4.23) 

and the Laurent series of these two operators coincide at cr = 0, since the contribution of the pole at = is, 
in both cases, given by the matrix Mq. Hence, we can define a meromorphic exponential dichotomy on M via 

{^%{x,y,(j) x>y>0 
l>^+(x,0,a)$i(0,y,a) x>0>y (4.24) 
$i(a;,y, cr) > x > y 

for X > y, and an analogous expression for <i>"(a;, y, a) for x < y. Note that, if we fix a; > and let y — ^ 0, we 
obtain from the first two equations in (4.241 the two expressions 

, $i(x,0,cr) 
$s(a; cr) = { +\ , , J 

^l{x,0,a)^tiO,0,a)^^l{x,0,a)il-P^{0,<j)) = ^l{x,0,a), 



which coincide due to (4.23). This completes the meromorphic extension of the exponential dichotomies on M 
for cr e Be(0). 

It remains to obtain pointwise bounds for these operators. It suffices to consider the case x > y, a,s the case 
a; < y is completely analogous. Recall that h±{a) = h^{a) + hp (a), where h^{a) is analytic in a, while 



±1 



where m° denote the entries of Mq. 



We have 

<I>^(a;,y,a) = ^l{x,y,a) - ^l{x,0,a)ht{a)^l{0,y,<j) - ^l{x,0,a)h+{a)^l{0,y,a) x>y>0 

" V ' ' 

(a) (i) 

^%x,y,a) = ^l{x,0,a)<^'^iO,y,a)-^^x,0,a)ht{a)^liO,0,a)<^UO,y,a) x>0>y 

^ V ■' 

(b) 

- $;(x, 0, a)h+{aW+iO, 0, a)$l(0, y, a) 
" ^ ' 

(ii) 

<^%x,y,a) = $l(x,y,a) + $^(a;,0,a)/i-(a)$l(0,y,f7) + $'i(x,0,a)/i-(a)$l(0,y,a) > x > y. 



(c) 



(iii) 



We now use the bounds (4.6) for the exponential dichotomies on that we derived in Lemma 4.1 First, we 
estimate the terms given by (a)-(c). For case (a) with a; > y > 0, we find 

W'f'ix, y, a)|U(y„) = ||<i>;(a;, y, a) - ^'^x, 0, a)ht (a) <i>l{0, y, a)|U(y,„) 
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Case (b) for x > > y gives 



\\<^'{x,y,a)\\^Y^-, = \\^'4x,0,a)^t{0,y,<y)-^\{x,0,a)ht{a)<^l{0,0,a)^t{0,y,a)U^Y,^) 

Finally, case (c) with > x > y gives 

\\^'{x, y, a)\\L(Y^) = \\^l{x, y, a) + ^^x, 0, a)K{a)^l{0, y, 

Next, we consider the meromorphic terms (i)-(iii). For case (i) with x > y > 0, we obtain 
$;(x,0,a)/i+(a)$t(0,y,(T) 

= U^4x,o,a)J2m°jUM{'^m,m^,y,'^)-)x 

= -<i>;(x, 0, a) [P-(0, 0) - Pl%0, a) + P-(0, a)] V m^UMim^, (0), ■)x 

= -^\{x, 0, a)[P;^(0, a) + 0(a)] ^ m°,!7,(0)(vl/,(y, a), ■)x 



= :1'3(y,o-) 



-<i>;(x, 0, a)F;^(0, a) ^ mO.C/,(0)(vI/,(2;, a), •)x + <^\{x, 0, a)O(l) ^ m°.C/,(0)(vI/ .(j/, a), •) 



E 
E 



*;(a;,0,0)P;^(0,0)TO°C/,(0) + O(e-3''-) 



(^'j(y,a),-)x + rest 



m°.{/,(a;) + 0(e-3''-) 



(*j(y,CT),-)x + rest. 



Note that the terms denoted by the Landau symbol O are all analytic in a. Equation (4.6 1 shows that j{y, a) 
<^\{Q,y,aY^j{Q) satisfies the bound 



\^,{y,a)\<KY,^- 



y>0. 



(4.25) 



We remark that, as a consequence of [22^, Lemma 5.1], (0, y, a)'^ is the stable evolution of the adjoint equation, 
and *j(y, (t) = $!j_(0, y, (t)-^^'j(O) therefore satisfies the adjoint equation. 

The term (ii) for x >0 > y can be treated similar to (i), and we obtain 

1 



:m°.L/,(a;) + 0(e-3'7-) 



{^jiy,a),-)x +rest, 



where 



|vl/,(y,a)| <i^^e-''r„(-)s^, y<0. 

Finally, case (iii) for > x > y gives 

^l{x,0,a)h-{a)^l{0,y,a) =Y, 



(4.26) 



-TO?-C/,(2;)+0(e-3')l-l) 



where \E'j(y,cr) obeys the bound (4.261. We summarize our findings in the following lemma. 
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Lemma 4.4 There exists an e > so that ( ^ . _? 0[ ) has meromorphic exponential dichotomies ^^{x,y,a) and 
<i>"(a;, y, a), defined respectively for x > y and x < y, for a G ^£(0) such that 



(4.27) 



(the minus and plus signs are for and <I>", respectively), where ^j{y,a) and ^^{x,y,a) are analytic in a. 



Furthermore, ^'(j/,cr) obeys the bounds ^^2^ and {4-26), while we have the bounds 



\\^'ix,y,a)\\nY^)<K< 



y - - 



'^i,t('^)2:p-i'i7,(<^)s/ 



<y < X 
y <0 < X 
y < X <0 



(4.28) 



for ^^{x,y,a), with analogous bounds for $"(2:, y, ct). Derivatives of ^^{x^y,a) and '^j{y,a) with respect to 
X and y satisfy the same estimates multiplied by [|cr| + e^'''^'] and \\a\ + e^'''^'], respectively, as operators in 



4.3 Pointwise bounds for the resolvent kernel 



We now use the asymptotic formulas for the exponential dichotomy given in Lemma |4.4| to complete the proof of 
Theorem [5] Lemma 4.4 shows that the meromorphic extension of the exponential dichotomy can be written as 



^''{x, y, a) = Ei{x, y, a) + E2{x, y, cr) + ^''{x, y, a), 



s,u, 



(3.231 of G, as 



where the i?i_2 terms come from the meromorphic components of (4.27). Thus, we may write the representation 
G,{x,y,(T-s) ^ Pi I \Ei+E2 + ^'']{x,z,(j)B{z)Vi{z,y,a-s)dz (4.29) 



-Pi 



Ei+ E2 + 
E1 + E2 + 



(a;, z, a)B{z)Vi{z, y, a; s) dz 
{x, z, a)B{z)Vi{z, y, a; s) dz, 



where B is defined in (3.201. We shall prove that the terms that involve <i>^'" may all be included in the term 
R{x,y,a;s) given in Theorem [s] while the terms involving £^12 lead directly to the terms Ei 2{x,y,a]s) of 
Theorem [5] Throughout, we will often suppress any dependence on s for notational convenience. 

We now consider (4.291 and begin by estimating the integrals that involve Each of these integrals can be 

estimated in the same fashion, and we therefore give details only for the integral that contains in the case 

(4.30) 



where a; > > y. For x > Q > y, equation (4.281 and Lemma 3.5 imply that 



< 



<^\x,z,a)B{z)Vi.{z,y,cT) dz 

\^\x, z, a)\\L(Y,^)e-'^^'-y\e^\^-y\ \B{z)Vi{z, y, a)|y„ dz 

1/2 



< C 



\\^^{x,z,a)\\l^y^^e~^^\^-y\dz^ e^^\^-y\\B{z)V,{z,y,a)\i.^dz 



1/2 



< C 



< c 



■t^^-y^B{-)Vi{;y,a) 



.''in. ''in 
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Derivatives of this integral with respect to y can be estimated as in (4.351 below, using again (4.281 and 
Lemma 3.5 Derivatives with respect to {x,t,s) lead to the same estimates as above, since they can be ac- 
counted for by estimating the derivatives of the dichotomies in L{Ym, Y^_i) instead of L{Ym)- In summary, the 
estimates for the terms in the exponential dichotomy transfer to the resolvent kernel and are captured by 
the term R{x, y, a; s) in Theorem |5] 

Next, we consider the integrals in (4.291 that involve the terms Ej. We focus on Ei as the term E2 can be 
treated similarly. Equation (4.271 implies that the two integrals in (4.291 that involve Ei can be combined as 
follows: 



Pi 



Ei{x, z, a)B{z)Vi{z, y, cr; s) dz 



(t) + Pi 



Ei{x, z, a)B{z)Ve{z, y, a; s) dz 



= --Ua:{x,t) / ^m^j {'i'j{z,a),B{z)Ve{z,y,a;s))^ dz 

^ J —00 ■ 



Assume first that y < 0. Using ( |4.8[ ) and the definition of ^^(z, cr), we have 



it) 



(4.31) 



(4.32) 



for z < 0. Thus, we may write 

/oo 
^m?^.(v|/,(z,a),6(z)y,(z ,y,a;s))^ dz 
-OO 



(cj-in(2:,cr),6(z)V^(z,y,cr;s))^ dz+ {^j{z,(T),B{z)Vi{z,y,a-s))^ d 



(4.33) 



^e"^.»(^'^^m°^. ( f e-''-(-)(-s^) (c,- in(z, a), e(z)l/,(z, y, a; s))^ dz 

■ \J~CO 

+^ '^)^l3{z)Ve{z, y, a; s))^ dz^ 



where iniV''^) ^ ^ each fixed (j/, cr), so that (4.31) becomes 



-uAx,t)J2^-'^-^''^yi-.^{y,af. 



We now claim that there are positive constants C and rj so that 

|Zi,in(y,^)| < C, \dylun{y,<j)\ < C [kl +e-''i^i 



(4.34) 



for all y < and all a with Rccr > — e. Furthermore, the same bounds are then true for s-derivatives of ^i4n(y, cr), 
since these are equivalent to i-derivatives, which are taken care of by the regularity of functions in Y^- This 
then yields the expressions (4.11 in Theorem |5] and the bounds (1.14) in Theorem [3] since we can expand the 
analytic term li^in(y,a) in a and subsume the cr-dependent part into R(x,y,a; s). 



The first estimate in (4.34 1 follows from (4.33 1 and Lemma 3.5 as in (4.30 1 above, since |e '^in(°')(^ y)\ < Ce^^'' 



It remains to show that differentiation with respect to y leads to exponential decay with respect to y. We estimate 
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here only the most dangerous component in (4.331: Lemma 3.5 and (4.321 give 
dy f e-''.n(-)(^-^)(c,,i„(z,a),S(z)y,(z,y,a;s))^ dz 

e--r„(-)(--!^) {cj.^^{z,a),B{z)dyVi{z,y,a-s))^ dz + 0(|a|) 







(/° 




\J —OQ 





= o 

V-/ —OQ 

Odal+e-^l^/l) 



and renaming rj gives the second estimate in (4.341. 



(4.35) 



5 Estimates of the Green's distribution and linear stability 

With the bounds on the resolvent kernel G■^, determined in Theorem [5] it is now possible to obtain the desired 
pointwise bounds on the Green's distribution by a simplified version of the analysis of |14l 133] . 



5.1 Bounds for the Green's distribution 

In this section, we prove Theorem[3] Corollary |3 . 6| states that the Green's distribution Q.^, is given as the contour 
integral 



27ri 



(5.1) 



that involves the resolvent kernel G*. Applying Cauchy's integral theorem and Theorem [5] which, in particular, 
states that G* may be meromorphically extended into the region {a € C : Recr > — e}, we first observe that 



(5.1 1 may be modified to 



27ri 



e"*[G,(x,2/,fT;s)](t) da, 



where f is defined in Figure [4] for some small constant r > 0. The fact that the key relation ( 3.14[ ) persists 
under analytic extension implies that the integration along the top and bottom pieces, [/i — i/2, — e/2 — i/2] and 
[— e/2 + i/2, fi + i/2], cancel; see Figure [4] We therefore have 



where 

r 



e 1 e . 

17. 

2 2' 2 



U 



Q*{x,t;y,s) 



ir, r — ir 



27ri 



'[G*(a;,2/,cr;s)](t) da, 



U [r — ir,r + ir\ U 



r + ir, — - + ir 



U 



This is a representation on a contour that corresponds exactly to the low-frequency part of the contour used 
to begin the arguments of [33l §8] and [TH §7], and we may therefore move the contours for the individual 
meromorphic pieces Ei, E2 and R of G* exactly as in these references. Since our resolvent bounds for G* as well 
as the initial contour F are the same as for the low-frequency estimates in the time-independent case treated 
there, we obtain the same bounds for Q^, that were obtained in [TH for the entire Green's function in the 
time-independent case. Note that the new term E2 that arises in G* for the time-dependent case has exactly 
the same form as the term Ei, except for the factor dtu in place of dxU, so it may again be treated in exactly 
the same way as before. We omit the contour estimates and refer the reader instead to ^S, §8] and [IT, §7] for 
details. This establishes the estimates (1.10)-(1.14l of Theorem|3] 
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i/2 


f 




(i) 




ir 




1 — 

-e/2 




1— 


[ 




(ii) 







Figure 4: Plotted is the contour f (solid and dashed lines) that is used to calculate the Green's distribution 
from the resolvent kernel G,. The contributions from integrating along pieces (i) and (ii) cancel, yielding the 
contour r (solid line). 



Finally, we recall that the singular parts Qq, Qi, and Q2 of the Green's function have already been estimated in 
Theorem |4] This completes the proof of Theorem |3] ■ 

We remark that the contour estimates carried out in the present time-dependent case are in fact somewhat 
simpler than those carried out in [Ml [23] for the time- independent case. The reason is that difficulties associated 
with high frequencies have been subsumed into the iterative parametrix-type construction of the resolvent kernel 
and were dealt with in Lemma |3.1| — an illustration of conservation of difficulty. Though we will not use this, we 
remark that the bounds we obtain for G* are nonsingular as < ^ 0^ and therefore somewhat better than those 
that hold for the low-frequency part of the time-independent Green's function. We shall use this observation 
later in the proof of Proposition |6.5| 



5.2 Proof of Theorem [T] 

From the bounds of Theorem [3] we obtain the equivalence of spectral and linearized stability as follows. 

As in the time-independent case, the pointwise bounds of Theorem |3] yield linearized orbital stability in CiH^, 
for any s, by standard — ^ convolution bounds; see [33j §9] and [21 §8]. This gives the sufhciency of 
spectral stability. Necessity follows from the theory of effective spectral projections, and the spectral expansions 
of the pointwise Green's function in terms of the effective spectral projections for {x, y) restricted to a bounded 
domain, that was developed in [331 §9] and [T31 §8]. We omit the details as they are essentially the same as in 
the time-independent case of those papers. 

This concludes the proof of Theorem [l] ■ 



6 Nonlinear stability 

We establish nonlinear stability by a modified version of the arguments used in [3 [12] to prove stability of 
general-type shock waves in the time-independent case. First, we will separate out the motion of perturbations 
along spatial and temporal translates of the shock profile from the rest of the time evolution by seeking solutions 

u{x, t) in the form 

u{x,t) = u{x-q{t),t-T{t))+v{x,t), Cit) ■.^{q,T){t). 

The decomposition of the Green's distribution into the terms £j, which capture the movement of perturbations 
along the shock profile, and the remainder Q, which describes the evolution along the characteristics, allows us 



to seek {(^,v) as fixed points of an appropriate set of integral equations that are similar to equations (2.5 1 and 



(2.6 1. An iteration scheme is then set up to construct solutions to these equations. To prove temporal decay, the 
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anticipated spatio-temporal decay properties of (C, v) are encoded in an appropriate set of template functions, 
which are used to control the norms of iterates. 



6.1 Fixed-point iteration scheme 



We now introduce the fixed-point iteration scheme by which we shall simultaneously construct and estimate the 
solution of the perturbed shock problem. We will use the notation 



■{x, t) = -{u,^,ut){x - q,t-T) 



and 



c 



-{ux{x - q^,t - T^),ut{x - q^,t - T^^:)) ■ {q{t),T{t)) 
= -u^{x - q^,t- T^)q{t) - ut{x - q^,t- n)f{t). 

We fix C = and consider the linearization about m^* . In particular, 

T^{y,s,t) = {Tri,TT2){y,S,t) 

will depend on C* through the profile u^* about which we linearize, though we will suppress this dependence in 
our notation. For later use, we note the important fact 



7r(2/, s,s) = 0. 



(6.1) 



Our starting point, similar to [3^], is the observation that, if u solves (1.2 1, then the perturbation v defined 



via 



satisfies 



where 



u{x,t) = u'^'+^'-^\x,t) + v{x,t) 



vt -C^'v^ [Q^' (C, v% + ^ 1^ at) + [R^' (C, v% + S^' (C, Ct) , 



(6.2) 



5«*(C,C) 



Vxx - [fu{u^')v]x 

(/(S«*+«) + A(#*+f)i; - /(u«*+« + «)) - 0{\v\') 
/.(S«-(x,t)) - Uu'^'+<{x,t)))v = 0(e-''l-l|C||«|) 



du'' 

C.+C(t) dC 



C = 0(e-''l-l|C||CI). 



The above asymptotics hold as long as \v\ remains bounded, by Taylor's Theorem together with exponential 
decay of d^u'^ in space, which follows from Hypothesis (HI) Note that, in the above, we have suppressed the 
time-dependence of C = C(^) ™ the superscripts for notational convenience. 

We can now describe the iteration scheme for (C*,C, w). For given C"^^ and C"^^('); 

Vg^^{x) ■.^uq{x)~u^" \x,0) 



(6.3) 



and define to be the solution of 



v"(x,i) 



g'^-\x,t;y,0)v^-'{y)dy + 



g'^-\x, t; y, s)Si''\C-\C-'){y, s) dyds 



(6.4) 
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where 5" is the decaying part of the Green's distribution CJ" — +^2 of the Hnearized equation 

around u^* . Further, set 



C{t) 



— 00 
t poo 



/O ^-co 

ot nOQ 



(7r"-i(y, 0, t) - 7r"-i(y, 0, oo))v^-\y) dy 

(7r"-i(y, s, t) - 7r^-\y, s, cx,))S'^'^\C'\ C-')(.y, -s) dyds 
{n^-\y, s, t) - ^»-i(y, s, 00)) [q^"" (C""\ «") + R^^'" (C"~\ (y, dyds 

<-i(y,s,oo) [gC-^(C"-\^;»)+i?C-^(^n-i,„")l (y,s)dyds 



(6.5) 



/•OO /"OO 



00 /'OO 



and 



n-1 



— 00 

00 /'OO 



7r--\y, s, 00)5'^""' (C"~\ C"~')(2/, dyds 



(6.6) 



00 
00 /'OO 



(y,s,oo) gf""(C"-\«") + i?^""(C"-\^;") (y,s)dyds. 



where tt" ^ = (tt" ^,7r2 ^) is associated with ^ and 82 ^ of the linearized operator at tt^* \ Formally, we 
can then define an associated iteration map T by 

We now clarify the relation between fixed points of T and solutions of the viscous system of conservation laws. 
Lemma 6.1 Under (6.4)-(6.6l, the function u'^ := vf"- +^ + w" satisfies the equation 

JC{t)-C-\t)) (6.7) 



9C 



with initial data u''{-,0) = ua{-) + {u'^' '+'^" '(°)(-, 0) - u?" '(-,0)), w/iere 



dt7Tr\y,o,tK-\y)dy 



00 

i /"OO 



dtnr'iy, s, t)^'^"" (C"~\ C"~')(2/, s) dyds 



^-00 



(6.8) 



^-00 



and 



T-{t) 



dt^^-\y,Q,t)v^-\y)dy 



— 00 

t poo 



a^^r i(y, s, t)S^'-\C-\C-')(.y. s) dyds 



+ / dtdy7Tl^-\y,S,t) Q^r\^n^\yn^ + RC:-\^n^l^^n^ (y,s)dyds. 



(6.9) 



In particular, u" satisfies (1.2) with initial data ua if and only if 

that is, if and only if (C", C") is a fixed point of T , in which case also C(0) = C,{oo) = 0. 
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Proof. Equations (|6.8| and (6.9l follow by differentiation of (6.5), recalling property (6.1) and the related fact 



that J^TTy{y,t,t)[Q, R]{y) dy — for functions Q and R of the type considered here; see [311 §4.2.4] for more 



details. From (6.5), we find that C^{oo) = and 

/>oo 

n^^-\y,0,t)v^-\y)dy- / / n"-\y,s,t)Si'~\C-\C-'){y,s)dyds 

Jo J-oo 

K~\y,s,t) \Qi"~\c-\v") + Ri"'\c-\v")] iy,s)dyds. 



C"W-C"(o) 



— oo 

t nOO 



Setting ^ = oo in this equation and comparing with (6.6), we find that C^{0) — C^(oo) = if and only if 



From (6.4 1 and (6.51 we conclude that 



g--\x,t;y,0)v^-\y)dy+ / / ^"-^(x, t; (C"-\ C"-')(y, s) dj/ds 

-OO ^0 J —oo 

ft pOO 



g--\x, t- y, s) b^"" (C"-\ «") + i?^"" (C"-\ (y, s) dyds+^ 

and thus, by Duhamel's Principle, 



Setting u" = w" + m^* "^"^ , we then obtain (6.7) by a straightforward calculation comparing with (6.2 1, with 
the claimed initial data 



(•, 0) + - Mo(-) + - S'^- 0). 



Note that the right-hand side of the above equation is equal to uq if and only if ^(0) = 0: if C" = C" ^ is a 
fixed point, this is true if and only if C"(0) = or, equivalently, C" — C"^^- ■ 



Remark 6.2 In (6.7 1 -(6.9), the values of (w", at time T depend only on the values for < t < T, and not 



on future times. By (6.4), we have, evidently, 



(6.10) 



6.2 Auxiliary convolution and energy estimates 

In this section, we shall collect various estimates that we will need to construct fixed points of the iteration 
scheme explained above. We begin by presenting two corollaries of Theorem [3j 

Corollary 6.3 Under the assumptions of Theorem^ and in the notation of Theorem^ the following holds for 
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j — 1,2 and y < 0: 

ItTj (y, s, i) — 7rj(?/, oo, s)| < C errfn | — ^ M for some a > 

\ M^{t-s) J 



s)lVAf(t-s) 



errin 



Symmetric estimates are true for y > 0. 

Proof. This follows from a straightforward calculation using and (1.141; see [131113] and [3 Remark 7]. ■ 

In the next corollary, we state a useful property of tt = (tti, 712). 

Corollary 6.4 ([ZJ) Under the assumptions of Theorem^ and in the notation of Theorem^ 

7r(y , s, 00) ^" ^5* \ y,s)dy^ ida2 . (6.11) 

Proof. This follows from the fact that both and are, for any fixed C, stationary solutions of the linearized 
equations. Hence 

g{x,t;y,s)ui{y,s)dy = ui{x,t), / g{x,t;y,s)ul{y,s)dy = ul{x,t). 



Because Ux and Ut are linearly independent and, under the assumption of spectral stability, £1 and £2 represent 
the only nondecaying parts oi Q{x,t;y, s), we have 

/>oo poo 

lim / ni{y, s,t)u'i{y, s)dy = lim / 7r2(?/, s, t)u^(?/, s)rfy = 0, 



which leads to (6.11 1 



Next , we investigate the dependence of tt on . 

Proposition 6.5 (Parameter-dependent bounds) Under the assumptions of Theorem^ and in the nota- 
tion of Theorem^ there exists a constant C such that 

\d^,dtA - C\dtT:\ 

l%(7i'(y,s,i) - 7r(2/,s,oo))| - C|(7r(y, s,t) - 7r(y,s,cx)))| 

\d^,dyn\ - C\dyTT\ (6.12) 

\c>c,{'^y{y, s,t) - TTy{y,s, oo))\ C\{'Ky{y, s,t) --Ky{y,s, oo)\ 

|%(9t9y7r| ~ C\dtdyiT\ 
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and 

\d^,d'^g{x,t-y,s)\ ~ C\d'^g{x,t-y,s)\ (6.13) 

for < |a| < 1 and y < 0, and symmetrically for y > 0, where by ~ we mean that the left-hand side obeys the 
same bounds as given for the right-hand side in Theorem^and Corollary \6.S\ above. 



Proof. Evidently, d^^ — {dx + dy,dt-\-ds)- The bounds (6.12) for tti and 7r2 follow by direct calculation, together 
with the estimate 



To obtain the estimates (6.13), we first indicate how to estimate s and t derivatives of Q, which we have not 
discussed so far. From the decomposition Q ~ X]j=o + ^* ^ ^ together with the bounds on Qj stated in 
Lemma |3.1| it is sufficient to estimate the t- and s-derivatives of 

g,{x,t-y,s) ^ £ = ^ i e''^'-'\R{x,y,a-smda, (6.14) 

where we have written the integral in terms of the unshifted time-coordinates, with initial time t = s, to capture 
the dependence on s — everywhere else in the paper, we shifted to i = 0, yielding the factor e'^* in place of e'^^*"'*-'. 



In (6.141, we see that t- and s-derivatives yield a factor ±cr, where they fall on e'^*^*~'*\ and otherwise contribute 
a t- or s-derivative of R. The additional factor a yields a contribution to the inverse Laplace transform that is 
smaller by a factor (1-|- (i — s))^5 than our estimate for \Q^, —£\. Likewise, derivatives that fall on R are harmless, 
since we established pointwise for those in Theorem [s] Similarly, using again the results from ^ the effect of 
additional spatial derivatives on (6.14) is a factor of order \a\ -l-e"'''^' for y-derivatives, while cc-derivatives do not 



change the bounds we obtain. Thus, for either temporal or spatial derivatives, the issue reduces to the estimates 
on Qj already carried out. ■ 

We shall also make use of the following technical lemmas proved in [7]. 

Lemma 6.6 (Linear estimates I) Under the assumptions of Theorem^ there exists a constant C such that 



f \g{x,t;y,0)\{l^ 

J — OO 




-3/2 dy 


< 


C'(^gauss H~ ^inncr H~ ^outcr)('^i ^) 


/ |7rt(y,0,t)|(l-f 

'J — OO 




-3/2 dy 


< 


C(l + t)-3/2 


/ |7r(y,0,t)|(l-f 

J —OO 




-3/2 dy 


< 


c 


^(y,0,t) -^(y,0,oo)|(l4 




-3/2 dy 


< 


C(l-h<)-l/2 



for < t < OO, where Q and n — (tti, tt2) are defined in TheoremwA 



Proof. These estimates can be established exactly as in |7] Lemma 3], using Corollary 6.3 and the explicit 
bounds for ^. ■ 

Lemma 6.7 (Nonlinear estimates I) Let 

e(y, S) := (1 s)^/2s"l/2(6lgauss + 6'inner + OonterfiV, s) + {I + s)"^(6'gauss + 6'inner + 6'outer)(?/, s) . 
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Under the assumptions of Theorem^ there is a constant C such that 



\gy{x,t;y,s)\Q{y,s)dyds < C(6'gauss + ^'inncr + 6'outcr)(a;, 



J-oo 
t poo 



J-oo 

OO POO 



\T:yt{y,s,t)\eiy,s)dyds < C{l + t)-^ 

|7r,(y,s,oo)|e(2/,s)d2/ < C{l + ty'/^ 

nT POD 

/ / \7ry{y,s,t)-TTy{y,s,(x)\e{y,s)dyds < C{l + t)-^/^ 

Jo J-oo 



for < t < OO, where Q and n = (tti, 7r2) are defined as in TheoremVA 



Proof. The estimates can be established as in [TJ Lemma 4], again using Corollary 6.3 and the explicit bounds 

for g. I 

Lemma 6.8 (Nonlinear estimates II) Let 



^2(y,s) 



— p-v\v\ 



(1 + ^) 



-3/2 



Under the assumptions of Theorem^ there is a constant C such that 



\Gy{x,t;y,s)\<^i{y,s)dyds < 0(6 

gauss H~ ^inncr H~ ^outcr)(^i ^) 



J-oo 

t pOC 



\nyt{y,s,t)\<i>iiy,s)dyds < 0(1 + t)-' 

\Try{y,s,t) -Try{y,s,(x)\^i{y,s)dyds < C{l + t)-^^^ 

nOQ nOO 

/ / \TTy{y,s,(x)\<^i{y,s)dyds < C{l + t) 

Jt J-oo 



-1/2 



and 



t I'OO 

/ \G{x,t-y,s)\<^2{.V,s)dyds < C{9 

gauss ^~ ^inncr ^~ ^outcr)('^; ^) 

J-oo 

t pOO 

/ \7Tt{y,s,t)\<f2{y,s)dyds < + 

J-oo 
pt poo 

/ / \TT{y,s,t)~TT{y,s,(x,)\<i>2{y,s)dyds < C{l + t)-^/^, 

Jo J-oo 

k(2/,s,oo)|$2(2/,s)dyds < + 



OO />oo 



ft J-oo 

for < t < OO, where G and tt are defined as in Theorem^ 
Proof. This follow as in T, Lemma 5] and fJW, Lemmas 4.2-4.3]. 



We shall also use the following auxiliary energy estimate, adapted essentially unchanged from [TSl UHl UHl I5T| . 
Let M be a solution of 

fill 

(6.15) 



for some given function j{t), and define the function v via 



(6.16) 
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Lemma 6.9 ( |15L I18L I31p Under the hypotheses of Theorem^ let uq £ , and suppose that, for < t < T, 
the suprema of \C,\ and I7I and the norm of v, defined by ( |6.15 1 and (6.16 1, each remain bounded by a 
sufficiently small constant. There are then constants 61^2 > so that 



(6.17) 



for allO<t<T. 



Proof. This follows by parabolic energy estimates similar to, but in fact much simpler, than those used to treat 
the case of partially elliptic viscosity in [15l [TSJ |3T] . Specifically, we write the perturbation equation for v as 

r-l 



Vt 



fu{u'^'+'^^'^ +Tv)dTV 



du{x) 



dC 



c. 



hit) + at)). 



(6.18) 



Observing that d"{du/d()\^^ (x) = 0(e ''l^') is bounded in the norm for \a\ < 3, we take the inner product 



in X of X^i^o ^x''''^ against (6.18 1, integrate by parts and rearrange the resulting terms to arrive at the inequality 



dtMUt) < -o\\dtv\\i. + c [Mis + i7(i)p + mr 

where 9 > 0, for some sufficiently large C > 0, so long as ||i'||_ff3 remains bounded. Using the Sobolev interpolation 

\\v\\l.<C-'\\dtv\\l.+C\\v\\l, 

for C > sufficiently large, we obtain 

dtMUt) < -oMh:^ + c (\\v\\h + m? + ICWP 



from which (6.171 follows by Gronwall's inequality. 



Lemma 6.10 ([19]) Under the hypotheses of Theorem^ let Mq := |1(1 + |xp)'^/'*uo(x) j] jif3 < 00, and suppose 
that, for < t < T, the suprema of |C| and \j\ and the norm of v, determined by (6.151 and (6.I6I, each 



remain bounded by some constant C > 0. Then there exists some M = M{C) > such that, for all < t < T, 



{i + \x\y^M^,t) 



< Me 



Mt 



Mo+ I (|CP + l7r)(r)dr 




(6.19) 



Proof. This follows by standard Friedrichs symmetrizer estimates carried out in the weighted norm. Specif- 
ically, making the coordinate change v = {1 + \x\'^)~3/4w, we obtain from (6.18 1 the modified equation 

»i 



Wt 



/„(m'^*+'^(*) +rw)dru; 



du{x) . 



IC 



Ai{t) + C{t)) 



(6.20) 



plus lower-order commutator terms of order 0(|w| -f \w\\v\ + |wj;||u|), which are bounded by M(|u;| -I- \'Wx\) by 
assumption, and similarly in the equations for for j = 1, . . . , 3. Likewise, J fu dr and its derivatives up to 
order two remain uniformly bounded by Sobolev embedding and the assumed bound on ||7;||/f.3. Performing the 



same energy estimates on (6.20) as carried out on (6.18) in the proof of Lemma 6.9 we readily obtain the result 



by Gronwall's inequality (indeed, somewhat better, thanks to parabolic smoothing). We refer to [19, Lemma 5.2] 
for further details in the general partially parabolic case. ■ 



Remark 6.11 Using Sobolev embeddings and equation (6.15), we see that Lemma 6.10 immediately implies that. 



z/||(l H- |a:;p)'^/'*Wo(2;)||_ff3 < 00 and if\\v{-,t)\\}j3, \C{t)\ and \"){t)\ are uniformly bounded onO <t <T, then 

\{l + \x\Y^Mx,t)\ and \{l + \x\Y^\{x,t)\ 
are uniformly bounded on < t < T as well. 
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6.3 Proof of Theorem d 



We are almost ready to prove the main theorem. We will need the following definitions and lemmas. Define the 
norm 

Kb, \am + tY^'Ur + IC(t)(i + t)\Lr (6.21) 

and the Banach space Si {C : IR ^ : |C|ei < oo}. We also define 

1 



-'gauss I ^mncr 



+ £'outor)(a;,t) 



v{x, t) 



(6.22) 



and the Banach space B2 '■— {i 



|w|b2 < 00}. The next lemma gives local existence of the integral 



equations ([6J])-(|6!9]) for (w",C"). 

Lemma 6.12 {H^ local theory) Under the hypotheses of Theorem^ let 

Ml := \\vo{x)\\h^ + \C-'\b, + ICr'l < 00, 



where \ ■ je^ is defined in (6.21). There exists some T = T(Mi) > sufficiently small and C = C{Mi,T) > 
sufficiently large such that, on < t < T , there exists a unique solution 



of (6.7l-(6.9l that satisfies 

|h;"(t)||«3, |C"(t)|<CMi. 

Proof. Short-time existence, uniqueness, and stability follow by (unweighted) energy estimates in u" similar to 
(6.191 combined with more straightforward estimates on ^" carried out directly from integral equations (6.8I and 
(6.9 1, using a standard (bounded high norm, contractive low norm) contraction mapping argument like those 
described in [31, §4.2.4] and [321 Proof of Proposition 1.6 and Exercise 1.9]. We omit the details. ■ 



Remark 6.13 A crucial point is that equations (6.7)-(6.9l depend only on values of (v",(") on the range 



t S [0,T]; see Remark 6.2 



The next result allows us to extend solutions provided they stay bounded in an appropriate sense. 
Lemma 6.14 Under the hypotheses of Theorem^ assume that 

{v\c)eLr{H'^)xLr 



satisfy (6.7 1 -(6.9 1 onO <t <T, and define 



sup 

xeR, s<£[0,t] 



K{x,s)\ 



eauss "I "inner ~r "outer 



){x,s) 



\cmi+s) 



(6.23) 



If/3{T), IIiJq '^llffs, and IC" ^je^ are bounded by some sufficiently small f3o > , then, for some e > 0, the solutior 
(w", C"), and thus (3, extends to < t < T + e, and (3 is bounded and continuous on < t < T + e. 



Proof. By (6.10 1, we have 



= ll«"(-,o)|U3, 

and Lemma 6.9 and smallness of both /3{T) and therefore imply boundedness (and smallness) of 

||u"(f)||_f/3 and \C{t)\ on < t < T. By Lemma [02] applied with initial data (w"(T), C"(r)), this implies 



existence and boundedness of in and in on < i < r + e for some e > 0, and thus, by Remark 
boundedness and continuity of/3onO<t<r-|-e. 



6.11 
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Lemma 6.15 For M > sufficiently large, AIq :— ||(1 + | • p)^/'*(uo(-) — u{-,0))\\h3 sufficiently small, and 



-^^IC" ^1 ^ 2CMq for some constant C > 0, there exists a Ci > sufficiently large so that solutions 



(w", C\C) of ^^-^^ exist for all t > and satisfy 
and 

|i'"|e. + IC"lei+M|C|<2CMo. 



(6.24) 
(6.25) 



Proof. Define (3 as in (6.231. We claim that, if we can show that 



(6.26) 



for some fixed C and > 0, for all time t such that the solution (w",C") e Lf{Hl) x C^" of ([6J|-([6!9| exists 
and 

m < ^CMo, (6.27) 



then we can conclude that the solution (w", C") in fact exists and satisfies (6.271 for all i > 0, provided 

2 

Mn < 



5C43C + 4) 



is sufficiently small. To see this, first note that, by (6.101, (6.3), and the exponential decay that follows from 
Hypothesis [(HI)] we have 



= \\U(i-lf' 



< \Wo-u\\h3 + \\u' 

< Mo + ci|Cri| 

f. 2ciC 



(6.28) 



which is small when Mq is sufficiently small. Defining T to be the maximum time up to which a solution (v", C") 
exists and f]{t) < /3q is sufficiently small, we can use the assumed bounds on ^"^i and (^"^i to apply Lemma 6.14 



and conclude that (u", C") exists up to time T + e for some e > and that /3 remains bounded and continuous up 



to T + e as well. Observing that (6.261 together with Mq < 2/(5C*(3C + 4)) implies that /3(i) < ^CMq whenever 



f3{t) < |CMo, we find by continuity that (3{t) < |CMo up to i = T + e as claimed. We may now repeat this 



process infinitely many times to conclude that the solution exists for all t > and satisfies (6.271 



We now show how (6.271 can be used to prove the lemma. By the definition of /?, (6.271 implies 



This proves part of equation (6.251. To complete it, we must show that 



i™i<^(i+r^/^ 



and 



IC 



< 



CMo 
4M 



(6.29) 



(6.30) 



Establishing (6.29) and (6.30) also proves that the integral equations for C,"^ and C" converge, and so we would 
also obtain, by Lemma 6.1 and the fact that («",(■") satisfies (6.7l-(6.9l for all t > 0, that {v^,(^'^,C") satisfies 
(6.4l~(6.6l as claimed. Finally, recalling (6.7) and applying Lemma 6.9 with 7 ^" — (^"^^, we obtain ( |6.24 l so 
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long as (6.261 remains valid, controlling ||u"||/f3 by integrating the right-hand side of (6.17) and using (6.271, the 
definition of (3, and the assumed bounds on C"^^. We shall carry out this last calculation in detail in equation 
(6.311, in the course of proving (6.261. Thus, it remains to prove (6.261, (6.29), and (6.30). 



We now establish (6.26) using (6.27). By Lemma 6.9 we have 



i/3 



< 



< c||i;"(0)||^3e-'" + c 



(K|i.+max{|C"|MC"-1^})(T)dr 



< C2(A/o'(l + 2ciC/Af)2 + (3CAfo/2)')(l + t)-^/^ 

< (ClAfo)'(l + ^)-^/^ 



(6.31) 



for Ci > sufficiently large and Mq sufficiently small, by (6.28), (6.27), and the definition of (3. With (6.2), 
(6.24), the assumption that \C\bi < 2CMo, and the definitions of /3 and | • Isi, we obtain readily 

IqC- + ii;C. I < c(/32 _^ 4C^M^){e + $i) (6.32) 

and 



(6.33) 



where Q and ^'1^2 are defined in Lemmas |6.7|]6.8[ Applying Lemmas 6.6 6.8 to (6.4), (6.8), and ( |6.9| , we 
thus obtain (6.26) as claimed. Likewise, we obtain (6.29) from (6.5 1 and (6.28), using Lemmas 6.6 6.8 and the 
definitions of /? and | • . 



It remains only to establish (6.30). This is more delicate due to the appearance of M in the denominator of the 
right-hand side and depends on the key fact that the estimate C" of the asymptotic shock location is to linear 
order insensitive to the initial guess To see this, decompose the expression (6.6) for C" into its linear part 



^"-i(2/,0,oo)<-i(2/)dy 
7i"lc.=o(y,0,cx))(u-?2^" ^){y)dy 



(6.34) 



7i'lc.=o(y, 0, oo)(wo - u){y) dy 



■■h 



n-l 



7i'lc.=o)(y,0,oo)wJ ^{y)dy 



=:-fo 



and its nonlinear part 



00 /"OO 



00 



^"-i(y, 0, ^)S<''\C-\C-'){y, s) dyds 
7r--\y,0,^) \Q^'~\u-,K) + Ri'~\C-\u-)\ iy,s)dyds 

L Jy 

By estimates like the previous ones, we readily obtain 

\II\ < 2c(2CMo)^ 



(6.35) 



which is < CM0/I6M for Mq sufficiently small. Likewise, \Ic\ < c|C*|Mo, by ( |6.28[ ), ( |6.12| , and the Mean Value 
Theorem, hence is ^ CMq/IGM for Mq sufficiently small (recall that we assume < 2CMq), and 

141 < c||uo - u\\li < C2||(l + \x\^f/\uo - u)\\h^ < C2M0, 
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hence is <C CMq/I&M for M > Q sufficiently large. Finally, Taylor expanding, and recalling | (HI ) | and (6.11 1, we 
obtain 



T ATI— 1 An— 1 

-'a — l,* 



7r|c.=o(2/7 0,oo) 



dy + 0(|C*n = 0(|C*P), 



which is also ^ CA/0/I6A/ for Mq sufficiently small (recall that we assume |C*| < 2CMo). Summing these terms 



up, we obtain (6.30) for Afg sufficiently small and M > sufficiently large, as claimed. This completes the proof. 



Proof of Theorem I 



We define 

\{vX)U \v\B^ + M\Cl {v,C)eBixR. 



Lemma 6.15 implies that, for r > small, M > sufficiently large, and 

A/o = i|(l + |-|T/^(«o(-)-^(-,0))|U3 

sufficiently small, the mapping T — T(C, (^^) is well defined from 

B{0,r) C 61 X X M. 

To establish the theorem, therefore, it suffices to establish that T is a contraction on 5(0, r) in the norm 
Indeed, we can then apply Banach's fixed-point theorem to find that 

r(c,c*) = (C,C*) 



has a unique solution, and Lemma 6.1 shows that the associated function u satisfies (1.2 1 with initial data Uq 



The stated decay estimates follow from (6.241, (6.25), and the definition of the norms in the spaces Bi and B2 



in (6.21) and (6.22) 



To show that T is a contraction, we need to establish the Lipschitz bound 

|r(c, c*) - T(c, C*)U < a|(c, C*) - (C, C*)|. 



(6.36) 



on B{0,r) for some a < I and some sufficiently small r. Assume that (w",C":C") satisfies (6.4)-(6.6) associated 
with iC^\C*^^), while (?)",C",C") satisfies (|6l])-(|6^ with (C""\C""^) replaced by (C"~\C"~^)- For each n, 
we define the variations 



Av 



"—{)"- u", AC" := C" - C 



n /-n /-n 



AC := C - C, 



and, likewise, define A^J" ^, Att" ^, AS"" ^, AQ" ^,andAi?" ^ in the obvious way. Using equations (6.4 1-( 6.6 1 
and the fact that - = A/".g" - /"Ag", we obtain 



Ai;"(a;,i) 



00 pOO ^ 



00 



Ag-~\x,t;y,Oyr\y)dy+ / ^"-^(x, t; y, 0)A<-^(y) dy 



00 



t /'OO 



Ag-~\x, t; y, s)S^''\C^\C-'){y, s) dyds 

^-00 

t /'OO ^ 

g^~\x,t-y,s)AS^-\y,s)dyds 

J -00 

t pOO 



J -00 

t nOO ^ 



Ag^~\x,t-y,s) [q^"'\v'-) + R^'~\C-\v'')\ {y,s)dyds 
/ / gi^-\x,t;y,s)[AQ"-^+AR-~^]{y,s)dyds, 

Jo J -00 
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Ant) = - 



oo 



oo 

t nOO 



dt7Tr\y,o,t)Av^"\y)dy 

> 

/ / dtA7Tr\y,s,t)si"-\c-\C-'){y,s)dyds 

Jo J-oo 

/ / dtnr\y,s,t)AS"^-\y,s)dyds 

Jo J-oo 

nt nOO 

/ / dtdyA7Tr\y,s,t) Q^"'\v'') + R^'~\c-\vn {y,s)dyds, 

Jo J-oo ^ ^ 

/ / dtdy7rr\y,s,t)[AQ"~'+AR^-']{y,s)dyds, 

Jo J-oo 



cso pOO 

Adt7T^-\y,0,tK-\y)dy- / ^t^r^y, 0, t)A<-i(y) dy 

-OO <J —oo 

dtA7T^-\y,s,t)S'^''\C-\C-'){y,s)dyds 

^-oo 

t pOO 

n — l/^. „ 4-\ A C71 — 1/ 



dt7rl^'\y,s,t)AS''-\y,s)dyds 



'0 J~ao 
i^t poo 



'0 J-oo 
pt poo 



10 J-oo 

and similarly for AC", where 



Now define 



sup 



9C* c.=cr- 

OdACr^e-"!^!!) 
|Aw"(j;,s)| 



^Acr' + odAcrTe-"!"!) 



+ |AC"(s)(l + ,s)| 



(6.37) 



dtdyAn^-\y,s,t) [Q^'-\v") + Ri'-\c-\v")\ {y,s)dyds 

pt pOO 

/ / dtdyn^-\y,s,t)[AQ"-'+AR"-^]{y,s)dyds, 

Jo J~oo 

/oo pOC 
A^"-i(y,0,oo)<-i(y)dy- / ^"-^(y, 0, oo)Az;„"-i(y) dy (6.38) 
-OO J —OO 

noo 
A7r"-i(y, s, ^)S'^^'\C-\C-'){y, s) dyds 
-oo 

noo 
7T"-i(y,s,^)A5"-i(y,s)dyds 
-OO 

noo 
A7r^-i(y,s,oo) (Q«"" (t;") + (C"-\ i;") (y,s)dyds 

-oo '- 

noo 
7r^-\y,s,^) [AQ"-i + Ai?"-i] (y,s)dyds, 
-oo 



(6.39) 



xeR, sG[0,t] \ (^gauss -|- ^innor ~^ ^outcr)(2^, s) 

Let f := |(AC"~\ AC-i)|* = |AC""^|i3i -I- Af|AC-i| be sufficiently small. From ( [OQ] ) and smallness of r and 
f we obtain immediately 

|AQ"-i + Ai?"-i| < C(rC(t) rP){e + $i). (6.40) 



Also, ( |6.32p and ( |6.33| hold with c(/32 + iC^M^) replaced by Cr^. We use ( |6.f2[ ) and ( |6.f3| of Proposition 
to obtain 

Att""^ ~ ttAC, < fTT, (6.4f) 
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and similar appropriate bounds for A^" ^ and its derivatives (of course, tt in (6.41 1 is defined at a point between 
C*^^ and C"^"'^). Next, using Lemmas 



6.6 



6.8 



in a procedure parallel to the one used in the proof of Lemma 



6.15 



we obtain 



Since ^(t) is finite and continuous in t, this estimate implies that 

Cr 



m< 



1 - Cr 



for a constant C that is independent of r and f. Now, replacing ^ in (6.40) with this bound, we substitute the 
result back into (6.371 and into the similar formula for A^". Notice that, with the exception of the first two 
terms, the terms in (6.371 are all quadratic in their source term, so giving us small enough bounds. Hence, using 
once again Lemmas |6.6H6.8] we obtain 

I AC" I < (^CMoP+^f + Crf^ {l + t)-\ 



(6.42) 



of which the two first terms in the right-hand side come from the first two terms of (6.371. Similarly, we obtain 

(6.43) 



|AC"| < (^CMo +^ + Cr^ HI + ty'^- 

We notice that {CMq + ^ + Cr) can be made arbitrarily small, provided that Mq and r are small enough and 
M is large enough. Next, we use (6.381 to bound A^", using basically the same method used in (6.34|-(6.35l, 
and therefore obtaining 

M|AC| < {CMo + Cr)r. 

This, together with (6.42) and (6.43), gives us (6.36) with a < 1, finishing the proof of the Theorem [2] ■ 



7 Summary and open problems 

In this paper, we considered time-periodic viscous Lax shocks u{x, t), which converge to constant time-independ- 
eni|^ rest states u± a,s x ~* ±oo. We showed that spectral stability of a time-periodic shock profile implies 
its nonlinear stability with respect to small initial perturbations that are smooth and sufficiently localized in 
space. Specifically, we proved that the corresponding solution converges with algebraic rate in L°° to an 

appropriate space- and time-translate of the shock profile u. 

The nonlinear stability proof followed the same strategy as in the case of stationary viscous shocks. First, the 
resolvent kernel of the linearization about the shock is extended meromorphically across the imaginary axis, 
and pointwise bounds are derived for this extension. These bounds together with analyticity make it possible to 
derive pointwise estimates for the resulting Green's function, using the inverse Laplace transform together with a 
careful deformation of the integration contour. Finally, a nonlinear iteration scheme that utilizes the anticipated 
spatio-temporal decay of perturbations closes the argument. The key difference from the case of stationary shocks 
is the time-periodicity of the underlying operators. To resolve this issue, we use spatial dynamics and exponential 
dichotomies to construct the analogue of the resolvent kernel and its meromorphic extension. One novel feature of 
our analysis is that we do not use Evans functions: instead, the geometric properties that ultimately determine the 
order of the pole of the meromorphic extension at the origin are encoded through Lyapunov-Schmidt reduction. 
Another advantage of our approach to meromorphic extensions through exponential dichotomies is that it is 
abstract and coordinate- free. This has several useful implications, which we shall now discuss. 

^Viscous conservation laws, ut + f{u)x = Uxx, do not support genuinely time-periodic homogeneous rest states. Thus, time- 
periodic shock profiles can only admit timc-indcpcndent asymptotic rest states 
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First, as already discussed in J|T| our results are also true for undercompressive, overcompressive, and mixed- 
type shock profiles, provided the notion of spectral stability is appropriately adapted and interpreted. Indeed, 
Lemma [4.2| essentially reduces the meromorphic extension to finite-dimensional problems that are analogous to 
the stationary case, and therefore well understood. Similarly, the nonlinear iteration scheme relies primarily on 
the template functions that describe the anticipated temporal decay and are identical to those of the stationary 
case. 

The fact that the Evans function is not used in our analysis opens up the possibility of extending our nonlinear- 
stability analysis to semi-discretizations of viscous shocks. Two-sided spatial finite difference approximations of 
viscous conservation laws lead to lattice differential equations. Their travelling waves satisfy functional differential 
equations of mixed type that are ill-posed as initial-value problems. More importantly, their stable and unstable 
eigenspaces are infinite-dimensional, and it is therefore unclear how Evans functions could be defined (we refer 
to |12l 116] for situations where this can be done through Galerkin approximations). It was, however, shown in 
[SI I13j that these equations admit exponential dichotomies. Thus, it is feasible that our approach could be used 
to establish nonlinear stability of semi-discretized shock profiles, thereby extending the analysis carried out in 
[1] for one-sided finite differences. 

We end this paper with three open problems. First, though we attempted to prepare the ground for a future 
analysis that addresses quasilinear, partially parabolic systems, we did not actually consider real viscosity here. 
We expect, however, that our main strategy should be applicable to such systems as in the related analyses of 
[IHJIIH]- Another open problem is to find concrete examples of time-periodic Lax shocks. While previous analyses 
have shown that time-periodic Lax shocks bifurcate at Hopf bifurcations from stationary shocks, no examples 
are yet known where these bifurcations occur. The last issue is of a technical nature. Our construction of the 
resolvent kernel relied on an iterative Birman-Schwinger or parametrix-type argument, which is quite involved 
and does not immediately give expansions or bounds of the resolvent kernel. As discussed in f|3] it should be 
possible to use exponential dichotomies directly to construct the resolvent kernel, but we have so far not been 
able to make this argument rigorous. 
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